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PREFACE 


In many Colleges of Engineering, the need is felt for a text- 
book on Differential Eyuations, limited in scope yet comprehen- 
sive enough to furnish the student of engineering with sufficient 
information to enable him to deal intelligently with any differen- 
tial equation which he is likely to encounter. To meet this need 
is the object of this book. 

Throughout the book, I have endeavored to confine myself 
strictly to those principles which are of interest to the student of 
engineering. In the selection of problems, the aim was con- 
stantly before me to choose only those that illustrate differential 
equations or mathematical principles which the engineer may 
meet in the practice of his profession. 

I have consulted freely the Treatises on Differential Equations 
of Boole, Forsyth, Johnson, and Murray. I am indebted to two 
of my colleagues, Professors N. C. Riggs and C. W. Leigh, for 
reading parts of the manuscript and verifying many of the 


answers to problems. 
D. F. CAMPBELL. 
CurcaGo, Ixt., 

September, 1906. 


PREFACE TO ENLARGED EDITION 


This book as it first appeared consisted of the first eight 
chapters as here given. The kindly criticism by a number of 
those teachers for whose use it was intended on the need of a 
discussion of equations, that occur in investigations in Mathe- 
matical Physics, other than those given in these chapters has 
induced me to add Chapter IX to the book. 


Af 


vi PREFACE 


In the preparation of Chapter IX., I have drawn freely from 
Professor Byerly’s Treatise on Fourier’s Series and Spherical 
Harmonics, from Professor Bocher’s pamphlet entitled Regular 
Points of Linear Differential Equations of the Second Order and 
from notes kindly loaned me by Professor Snyder of Cornell 
University. I have also consulted Heffter’s Treatise on Linear 
Differential Equations with one Independent Variable. 

_ To those teachers who have sent me their criticism of the 
book in its original form, as well as to others who have cordially 
received it, I am under the deepest obligations. 


D. F. CAMPBELL. 
Cuicago, ILx., 


June, 1907. 
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A ee COURSE 


DIFFEREN TAL EQUATIO 


CHAPTER I 
INTRODUCTION 


a oblo, We shall here adopt the following : 

If to every value of x there corresponds one or more value 
f(«), then f(«) is said to be a function of z. 

This definition includes a constant as a function of 2, 
f(z) is constant, then for every value of x, f(x) has a vs 
namely, this constant. 

A definition of a function of two variables is the follow 

If to every pair of values of two variables # and y the 
responds one or more values of f(«, y), then f(a, y) is said 1 
a function of # and y. 

This includes a constant or a function of one variable 
function of x and y. 

A function f(2) of one variable « is single valued when for 
every value of « there is one and only one corresponding — Vv 
of f(z). | 

A function f(w) of one variable x is continuous for a 
x = aif f(a) is finite, and 


Hal SiG xi h)| ue ee [ f(a ae h)| = f(a). 


A function f(2, y) of two independent variables x and y is 
single valued when for every set of values for w and y there is 
one and only one corresponding value of f(z, y). 

: 1 
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A function f(x, y) of two independent variables x and y is 
continuous for a set of values c= a, y = b if f(a, b) is finite, 


and 
limit 


w=] sat h, b+ ‘| = f(a, 6) 


no matter how h and & approach zero. 

The following definitions are given in almost any work in 
calculus : 

If f(x) is a single valued and continuous function of 2, given 
by the equation y = f(x), then 

Aw and Ay denote the increments of x and y respectively, 


dy li mit Ay 
dx 4%=0} Ax)’ 


If f(z) is single valued and continuous, and dy/dz is contin- 


uous, then 
Ey eee 
dx? ~ dx \ dx }* 


In general, if f(#) is single valued and continuous, and the 
preceding derivatives are all continuous, then 


dy d dy 
dx* ~ dx \ da* }* 
If f(a, y) is a single valued and continuous function of two 
independent variables x and y, given by the equation z = f(a, y), 
then 0z/Ox is the derivative of z with respect to x when y is held 


constant; Oz/Oy is the derivative of z with respect to y when is 
held constant. 


2. Inasingle valued and continuous function f(a) of one vari- 
able a, given by the equation y = f(x), whether 2 is the inde- 


pendent variable or a function of some other variable or Nala 
we have 
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dx = d(dx); dx = d(d’x); ---; dw = d(d™"2); 

dy = d(dy); dy = d(d’y); ---; d°y = d(d"'y). 
Definitions. The differentials dz, d’x, d’x, +--+, d’x, or 


dy, d’y, d’y, ---, d"y are called the first, second, third, ---, nth, 
differentials respectively. 


3. Derivation of d’y and d*y when no assumption is made re- 
garding x being independent or a function of some variable or 


variables. 
d'y = d(dy) =d ( iar ) = = eae +o ads 


2 dy 1. , Yn 


dy Y and'x 4. 
= 5 da! + 874 dud'a Pa. d*x. 

By taking differentials in succession any differential may ulti- 
mately be found. 


4, In the differentials of the preceding article, if x is an inde- 
pendent variable, it can be assumed without loss of generality, 
that Az, or what is the same in this case, dz, is constant. That 
is, it can be assumed that x changes by equal increments. Under 
this supposition, therefore, d’x and all higher differentials of « 
can be taken zero. Therefore, under this supposition, 


2, — 4Y 
dy = 73 —~ da’, 
By = Shae. 


The place which a derivative or differential occupies in the 
succession of derivatives or differentials indicates the order of the 
derivative or differential. Thus, a second derivative or differ- 
4 ential cme to be of the second order, a third of the third order, 
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5. The only functions usually considered in elementary works 
in calculus are functions of a real variable. Such functions with 
one exception are the only ones considered in the following pages. 
The exception is ¢& where z is a complex quantity. 

The student is already familiar with the definition of e* where 
xis real. He is, however, probably not familiar with the defini- 
tion of e when z is a complex quantity. A definition of this 
function will now be given. 


The infinite series 
2 23 


z 
Pgs ea 


where z is a complex quantity, can be shown to have a determi- 
nate, finite value for every value of z. It also reduces to the 


infinite series 
3 


1 tong ces . 
ae EP 8 


when z becomes real and equal to x, and this series, it will be re- 
membered, is equal to e* for all values of x. It therefore appears 
that the infinite series in z would be satisfactory as a definition 
of &. We shall define e* by saying that it is equal to the infinite 


series 
3 


a ae 
pad ig Pig ae 


for all values of z. 
From this definition, the following theorem can be established : 
Theorem. Ifz—x+y where x and y are real, andj = ¥—1, 
then 


e = &(cosy+)siny). 
Proof. 


& 
2g? 


3 
‘eal +8 iy tig +7 by definition, 


=14 (+a) 4 SGU, Cow... 
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Consider all the terms containing 2”. These are found from the 
terms 
(+ y)" i (e+ yy)" 
e eat 


— oe 
They are 


or 


ye 3 4 be 
lay — ~245+8_-.]. 
al 2 foe asa 


Separate the real terms from the imaginary and there results, 


Hh ee 
[r [cos y +‘j sin y]. 


or 


Let rv take all positive integral values in succession from 0. 
In this way we get all the terms of the development ¢. Then 


= (Pos Se 25 oe [cosy + j sin y] 
27 8 i 


= (cosy +Jjsiny). 
The theorem is therefore proved. 


Exampies. ¢ 7” — e“*(cos 3% +) sin 32) 


EXERCISES 

1. Given y = log a, find dy, d’y, d’y: 

(a), on the assumption that x is the independent variable ; 

(6), making no assumption with regard to a. 

In the results of (b), substitute x = cos @ and show that the 
results are the same as those obtained by first substituting the 
value of w in log # and then taking the differentials. 

2. Given y = & where x = cos 6, express dy, d’y, d*y in terms 
of 6 without substituting the value of « in the equation y = e*. 
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3. Given y = log where « = sin 6, express dy, d’y, d’y in 
terms of 6 without substituting the value of 2 in the equation 
y = log x. 

4, Prove that &™ = e*(cosy —jsiny). 

5. Prove that ete7*# — e*t**o", where x, y, z and w are 
real. 


ANSWERS 
Ci aes ae Qda3 
1(a). dy =—-3 fy = — ; y= eo 
ada — da? vda— 32dad’?x+-2dx* 
(6). dy ==; Py= ae ke Dy= 2 re 


2. dy = — sind 6° dé ;sx 
d’y = — sind e°°d’é 4. (sin’6 — cos0) Ee °d@ ; 
d’y = — sin e°*d?6 + 3(sin’@ — cos6) es *déd?6 
+ sin6(1 + 3 cosé — sin’6) ccs ?d@. 
3. dy = cot@ dé; d’y = coté d’6 — cosec’d dé ; 
d’y = cot d°6 — 3 cosec’é déd’6 + 2 cosec’é coté dé. 
6. Definition of differential equation. A differential equation 
is an equation involving derivatives or differentials with or with- 


out the variables from which these derivatives or differentials are 
derived. 


The following are examples of differential equations : 


dy\? 

(33) =a (1) 
dy | (dy\’ 
wat (GL) ty=0. (2) 
a dy, 
oat ya t+ yf = 0. (3) 


dy\? dy dy dy (dy 
3(53) - dx da? da? z) =o C2 
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Oz Oz 0 
Bast By (5) 
Ovz O*z ar 
dat t By oie oy 


7. In examples (1) to (4) inclusive of the teens article 
it will be noticed that differentials enter the equation only in de- 
rivatives. It is conceivable, however, that there might be an 
equation containing differentials other than those in the deriva- 
tives, as for example, 

= + 2x = + y=”, 
but there is no need of entering into a discussion of such equa- 
tions, and we shall not do so. In what follows, we shall assume 
that if the equation is written in differential form, the differen- 
tials can all be converted into derivatives by the process of 
division. 

8. Classes of differential equations. Differential equations 
are divided into two classes : ordinary and partial. 

An ordinary differential equation is one in which all the 
derivatives involved have reference to a single .independent 


variable. 


A partial differential equation is one which contains partial 
derivatives and therefore indicates the existence of two or more 
independent variables with respect to which these derivatives 
have been formed. 

Thus, in Art. 6, equations (1), (2), (3) and (4) are ordi- 
nary differential equations, and equations (5) and (6) are par- 
tial differential equations. 

Chapters I to VII inclusive are devoted to a discussion of ordi- 
nary differential equations. Chapter VIII contains a short 
treatment of some partial differential equations. 

9. Order and degree of a differential equation. The order of 
a differential equation is that of the highest derivative or differ- 
ential in the equation. 
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Thus, in Art. 6, equations (1) and (4) are of the third order, 
and (2) and (3) of the second order. 

The degree of a differential equation is the degree of the deriv- 
ative or differential of highest order in the equation after the 
equation is freed from radicals and fractions in its derivatives. 

Thus, in Art. 6, equation (1) is of the second degree, equa- 
tions (2), (3) and (4) of the first degree. 


Solutions of a differential equation. Let us consider the 
differential equation in each of the two following examples, and 
see if, from the equation, we can get a relation connecting x and 
y and not involving derivatives, such that, if the value of y in 
terms of x be substituted in the equation, the equation is satisfied. 


EXAMPLE 1. dy = 2. 
dx 
By integration, we get 
la 
y = 3 +6 
EXAMPLE 2. a +y=0. 


_ dy am ; 
cs a vVe—y 
Awa ° 


.yst vesin(e+e¢), or y=+-vecos Gare): 


In example 1, if $x° + ¢ be substituted for y in the equation, 
there results «7 = x’. The equation is therefore safisfied. 

In example 2, if + Vesin (w + ¢), or + Ve cos (# 4 ¢,) be 
substituted for y in the equation, there results, in the first case, 


INTRODUCTION 9 


= Vesin (+ ¢,) + Vesin (2 4 ¢,) = 0, and in the second case, 
= Vecos (# + ¢,) + Ve cos (x 4 c,) = 9. In either case the 
equation is satisfied. 

Definition. A solution of a differential equation is a relation 
between the variables of the equation and not involving deriva- 
tives, such that if the value of the dependent variable be substi- 
tuted in the equation, the equation is satisfied. 

Thus, y = 42° + ¢ of example 1, and y = + Vesin (x +. ¢,) 
of example 2, are solutions of the equations. 

In this book we shall not concern ourselves with the question 
of whether every differential equation has a solution but shall be 
content with finding solutions in the few special cases discussed 
here. 

11. A solution of an ordinary differential equation may be 
one of three kinds: general, particular and singular. 


A general solution is one which contains arbitrary .constants 
equal in number to the exponent of the order of the equation. 

Thus, in example 1, Art. 10, the number of arbitrary con- 
stants is one and the exponent of the order of the equation is 1, 
and in example 2 of the same article the number of arbitrary 
constants is two, and the exponent of the order of the equation 
is 2. In either case the solution is the general solution of the 
equation. 

A particular solution of a differential equation is a solution 
obtained from the general solution by giving one or more of the 
constants particular values. 


Thus 
y? y 1 ye 5 
Y= 3 Ur== 3 — or y= 3 —v, 
of example 1, Art. 10, or y = sin, y=2sinz, or y= —3 cosa, 
of example 2 of the same article, are particular solutions of the 
equations. 


A singular solution of a differential equation is a solution with- 
out arbitrary constants which cannot be derived from the genera} 
solution by giving the constants particular values. 

2 
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Singular solutions will not be considered in this book. 


12. A solution of a differential equation is not a general solu- 
tion unless the constants are in number equal to the exponent of 
the order of the equation, and cannot be reduced to a fewer 
number of equivalent constants. 

Thus, y = ce*t*, c and @ arbitrary constants, although it con- 
tains two arbitrary constants, is not the general solution of a dif- 
ferential equation of the second order, as can readily be shown. 
The solution y = ce*** is the same as y = c*e*. Now c* is equiv- 
alent to only one arbitrary constant because an arbitrary con- 
stant can have any value and thus all the particular solutions 
got by giving ¢ and a all possible values can be obtained. There- 
fore y = c%e* is equivalent to a solution y = Ae’, A arbitrary, and 
cannot therefore be the general solution of a differential equation 
of the second order. 


13. Let y =f,(%), y=f,(2), +++, y = f,(2) be solutions of a 
differential equation. 


Definition. If the c’s cannot be chosen, not all zero, such that 
f(z) +¢f,(@) +--:+e¢,f,() is identically zero, then the 
solutions are said to be linearly independent. 

Thus, y = + Vesin (x +. ¢,) and y= + Vc cos (# 4 ¢,) of ex- 
ample 2, Art. 10, are such that no values ¢, and ¢,, not both 
zero, can be chosen such that + ¢, Vesin (x+¢,) +e, Ve cos (#-+¢,) 
is identically zero. The solutions are therefore linearly inde- 
pendent. 


14. Derivation of an ordinary differential equation. Let 


(2, Y ¢,) = 0 (1) 


be an equation containing « and y, and the arbitrary constant ¢,. 
By differentiation of (1) there results 


Of Op dy 
Ca Cie (2) 


Equation (2) will in general contain ¢,. If between (1), and 
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(2), ¢, be eliminated, the result is a differential equation of the 
first order of which $(«, y, c,) = 0 is the general solution. 


EXxampLe. Find the differential equation of which 


is the general solution. f r\ 


ras 


Eliminate ¢, between the equations. Therefore 5 5 -m+t Z6,€ 


dy 4+ 2ay= me . Lou 


daz ad 
} 


is the differential equation of which 


m : 
f= 9 Si ee~ 
is the general solution. 

Sometimes the arbitrary constant is so involved that it disap- 
pears in the equation which results from the differentiation. In 
such a case this equation is the desired equation. 

Exampee. Find the differential equation of which y’ = 2¢,x 


is the general solution. 
Divide both sides of the equation by 2. 


1° 


2 
a Y | 2e 
xv 


By differentiation there results 


ae 


da a 


which is the desired differential equation. 


Let 
(2, Y, %, ¢,) =A (1) 


be an equation between 2 and y, and two arbitrary constants ¢, 


and ¢,. 
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By differentiation of (1) there results 
Op Oddy 0 (2) 


dat By de 
Equation (2) contains dy/dz and will in general contain ¢, and 
c, also. Eliminate one of the constants between the two equa- 
tions. Suppose the constant ¢, to be eliminated. The resulting 
equation contains dy/dw and in general 2, y ande,. Call it 


: dy 
W (a, Y, ae ¢,) = 0. 
By differentiation there results 


yey Sey (3) 


Ox oy dx 


Equation (3) contains d’y/dx’ and will in general contain ¢,. 
Eliminate c, between (2) and (3). The result is a differential 
equation of which $(, y, ¢, ¢,) = 0 is the general solution. 


Exampe. Find the differential equation of which 


C, 
Y = Oe + i 
is the general solution. 
Differentiate yY=eu+t 2 
x 
ee dy =c,— Co 
da ee fee 
Eliminate ¢,. 
dy  2¢, 
Differentiate. 
ek 
Oe eae 
Eliminate c, between 
dy Qe, dy  2¢, 
dat = and yea = 
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aot ot 
i4e de J =% 


which is the desired differential equation. 


15. It is seen from the preceding article that one constant can 
be removed after each differentiation. From this it would be 
expected that, starting with the differential equation, an arbi- 
trary constant might be introduced every time the order of the 
differential equation was lowered by unity. Then, since lower- 
ing the order of a differential equation of the mth order by unity 
nm times would result in a solution of the differential equation, it 
would be expected that a solution would contain not more than 
n arbitrary constants. 

It isa theorem that a differential equation cannot contain a 
solution having more arbitrary constants than the exponent of 
the order of the equation unless the constants are such that they 
can be reduced to a fewer number of equivalent constants. This 
will be assumed without further discussion. 

It is also a theorem that a differential equation cannot have 
more than one general solution. This theorem will be assumed 
without discussion. 


16. A general solution may have various forms but there is 
always a relation between the constants of one form and those of 
another. Thus, the general solution of example 2, Art. 10, may 
be written y = A sinw + B cos « instead of y= + Yesin (x+¢,). 
This latter form of solution is y= + ve cos ¢, sin w+ Ve sin ¢, cos 4, 
so that A = + ve cose, and B= + vesing,. 


EXERCISES 
1. Determine the order and degree of the following equations. 


ik d 2 A dy 2) 3 
(a). oi e{14 (51) } (c). sae 1+ (2) } 


ad? oa d Ou Ou 
(Db). fate aat a +y= 0. (d). GE = As 


14. SHORT COURSE ON DIFFERENTIAL EQUATIONS 


In each of the seven following exercises determine the differ- 
ential equation of which the given equation is the general solu- 
tion, given that ¢,, ¢, and ¢ are arbitrary constants. 


2. y = ¢, Sin mx + ©, Cos mx. BL ofan Cb Oe. 
3. v = ¢, cos (mt + ¢,). 6. vy = Ce 3 Oe 


aoe —¢,) (yy — 6) =m. 1. ft — Zee = 6" = 0. 


2) aa 
8. y+. 2°)? = ¢ + log, Tyee § 
9. Show that 


is a solution of 


10. Show that 
1 5 
4y = 35 $a + ae 


is a solution of 
dy 


d’y i 
2 ne — ss 
Ras ei Oy ag 


11. Show that 


y = pre + o67* + oe * 
is a solution of 


12. Show that 


c 

v= a ++ C, 
is a solution of 

dv Qdv 0 

det Tp dp 
ANSWERS 

d? dy 
2, cas — my. 3. aa = — mo. 
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dy : m d’y 
z 1+ (i) oy dy ae 
Lay 
dy dy  (dy\° d'y 
ea "det dx (2) : : dx? 
dy dy\? dy 
7. y— 20 —y (Zt) =0 we 


15 


dy 
on Bee 
uo L 
1 gtd = a t-ea*) 


CHAPTER II 
CHANGE OF VARIABLE 


17. Interchange of Variables. It is sometimes desirable to 
transform an expression involving derivatives of the function y, 
in y = f(x) where 2 is the independent variable, into an equiva- 
lent expression involving derivatives of the function 2, given by 
the same equation, where y is the independent variable. 

The formulas for such a transformation can be readily estab- 
lished as follows: 


dy limit | Ay] limit [ 1 J 1 (1) 
dz 4z=0} Ar | 4y=0] Ar | dx 
Ay | dy 
aye d / dy 
dx? ~ dx \ dx 
_# (ay) dy since @_ ty 
= dy\ dx) de’ ° da = dy’ de 
halal 1 
a, ce ee by substitution from (1). 
dy] dy 
dx ax 
dy’ 1 dy? 
i: Be D faene 2) 
dy} dy (o) 
dx 
Gy da [dy d dy? i 
Cine dx 
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a dx 3 d’a\* (da\2 dx dx 3 aay. 
_ May) ap T° \aep) \ty) af ty Tay 
ee dx 7 = dz 5 
dy ( dy 
The method of procedure for higher derivatives is evident. 


The transformations to which these formulas apply are called 
change of the independent variable or interchange of variables. 


~@ 


ExampLe. Change the independent variable from a to y in 


the equation 
3 By\* dy dy Py (dy\* 0 
OP ate dak = de ae 


Substitute from (1), (2) and (3). 


Gia y? dx dx d?z \? 

Bee nO) 
(de® | dz dz\° 
Pen eal) 

dx 


: ( ie d«\ dz 
"\ dy a ia) dy ~ 
18. Change of the dependent variable. Suppose that y is a 
function of x and at the same time is a function of some other 
variable z. The derivatives of y with respect to x can then be 
exptessed in terms of derivatives of 2 with respect to «. 
As a function of z, let y = $(z). Denote differentiation with 
respect to z by primes. Then 


dy dydz do(z)dz 
de dedz” dz dx 


d 
= ¢'(z) . 


dz? da\dx 
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ay d dy yr ” dz dz , 
fy =o (4) =0°' (GE 5)’ 4 88") Gg + 8 O) ae 


Similarly for higher derivatives. 
The above transformation is called change of the dependent 
variable. 


Exampte. In the equation 


2 Uy dy" 2 dy 
+9) F4— @y—1) (Gt) + de ty) Ge = 0 
change the dependent variable from y to z where y = tan z. 
dy 


dz 
aE a= Gay 


dx dx 


d’y dz dz 
dg = 2 See" etane (7) + sec as 


Substitute in the equation. 


2 


dz\e dz dz 
*, 2 sec* stano(7:) +sec* 27a — (2 tan z—1) sec +(Z) 


dz 
4 ——_—=— er 
+ 32 sec a7 = 0. 


dz dz dz 
dat (F ) Te rige tlt 

19. Change of the independent variable. Suppose that y is 
a function of x where x is a function of some other variable z. 
The derivatives of y with respect to 2 can then be expressed in 
terms of derivatives of y with respect to z. 

As a function of z let z = $(z). Denote differentiation with 
respect to z by primes. Then 
dy dydz dy1 dy 1 
dz” dzdx” dzdu~ dzd¢'(z) 


dz 
Py nd) = d (2 1 \z spose! d’y $’'(z) dy 
da* dx\dx}) ~~ dz\ dz $'(z)/du {¢'(z)}*dz?~ {#/(2)}* dz 


dz 
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dy d (33) Pa | 1 d’y oe) dy | 1 
er acd, ef da|_ {¢'(2)}* dz* {$'(2)}° dz | dx’ 
dz 

eg ee) ey 


¢’ z) 3 de? {'(z)}4 dz? 


_ # OH") — 6"@ Hedy 
{v@y} dl 


Similarly for higher derivatives. 


The above transformation is called change of the independent 
variable. 


ExampLe. In the equation 


d’y 2— Oy y 
eg mss re ams Baa a 


change the independent variable from x to z, where x = cos z. 


dy dy dz dy 
ieee cite Fad gs ae FS 


ay d (dy d dy 
dat ~ aes = de (- eg Le 


d 
= — cosec’z cot god + cosec? z 
2 


d. 


2, 


dz” 
Substitute in the equation. 


2, 


d d 
*. cosec” z a cosec” z cot z of + cosec’ z cot 2 e + cosec’z y = 0. 


When changing either the dependent or independent variable 
to a third variable, it is better to work out each derivative in the 
particular case considered rather than use the derivatives ex- 
pressed in the general case as formulas. 
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EXERCISES 
In each of the four following exercises, change the independent 
variable from x to y. 


d’y dy d’y dy \? dy \? 
dy \? | 
By ee) lee 
aga Uh en 
dx’ 


d’y dy d?y \? dy (dy\? dy\* _ (dy * 
ee 8 (Ze ee dix ae Ao hee di) = 
In each of the two following exercises change the dependent 
_ variable from y to z. 


5. +a) h— (FH) 20 4H vey, 


- where y = tan z. 


d’y dy d? 
2 Aare WTS 
Cae a (ay + Bay") S 


+ {2 (3) + ay + 32%y'} eH + xy’ = 0, 
where y = é’. 
In each of the four following exercises, change the independent 
variable from «x to z. 


is eT ateg ty = Te 
d 


dy 


2 
q 
8. (1 — a) 54 — woe = 


0, where 2 = sin z. 


9. Poa oa +v= 0, where wv = &. 


dy 2x dv v 


aD de +14 adet (4 a) 


5 = 0, where # = tanz. 
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11. Transform the formula for radius of curvature, 


+(e) | 


=o aa 


dx? 


into polar codrdinates, the equations of transformation being 
x =7rcos6, y=rsin O. 


ANSWERS 
dex fio Cn Ox 
dy de dx dix 
dy 
dz dz ae d*z d’z A ee 
Ds de ~ de = SD * 6. di ~ te gga t+ Bt Gg, += 0: 
a? P 
7. aa ty=0. 8. = 0. 
du dv dy 
9. Fa qa tr=?. 10. qa tv=2. 
dO 
iE. P= ak 


CHAPTER III 


ORDINARY DIFFERENTIAL EQUATIONS OF THE FIRST ORDER 
AND FIRST DEGREE 


20. An ordinary differential equation in one dependent vari- 
able, of the first order and first degree, may be represented by 
the equation 


Mdx + Ndy = 0 


where M and WN are functions of x and y and do not contain 
derivatives. 

The equation Mdz + Ndy = 0 cannot be integrated in the 
general form. There are certain particular forms of it, however, 
which can be integrated. Some of these will now be investigated. 


21. LINEAR DIFFERENTIAL EQUATIONS OF THE FIRST ORDER 


Definition. An ordinary linear differential equation of the 
first order is an equation in the form 


dy |p 
Peg eee ek! 


where P and @Q are functions of # and do not contain y or deriv- 
atives. 
The general solution of the equation 


dy 
tty = 


can be found as follows : 
Multiply both sides of the equation by eJ Pi, 


tn Ot es Pie On) 


22 
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If the substitution w — ye/?* be made, the left hand member of 
the equation reduces to du/dz. 


eure ee = [Qed do + ce S Pax (1) 


which is the general solution of the equation. 
In the original equation, if P is zero, the equation reduces to 
the familiar form dy/dz = Q, and the general solution is 


=e +.f Qde. 


If Q is zero, the equation becomes 


and the general solution is y = co SJ Paz, 
When Q, in the equation 


di 


dx 


is zero, the equation is called the ordinary linear differential 
equation of the first order with the right hand member zero. 


ExampLe. Find the general solution of the equation 
dy 1 2 
: Jee 
Multiply both sides of the equation by e7 2”. 
A dy S ite 1 Ss 2 S <4 
ome da + é y= ve - 
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1 
Let wu = ae 


1 
Ne sel e™, its [vel 2 "de + ©. 


1 
Now Be zo” — ge? _ x, 
a" 1 /2* et 
us fodepe= pte. syat(Fte)=F4+$ 


It is usual to solve an ordinary linear differential equation of 
the first order by substituting directly in formula (1). Thus, 
in the above example, formula (1) becomes 


1 
y= os aoe eS s*du + co Sa” 


oe 
de 
22. EQUATIONS REDUCIBLE TO THE LINEAR FORM 
A form easily reducible to the linear form is 


dy 
de t PY = Oy" 


where P and Q are functions of x and do not contain y or deriv- 
atives. 
Divide by y”. 
By Et. 8. 
dix 
Let yj"? = u. 


., dy, du 

*.(—n+1)y pain 
1 du 

See ie Q 


: = sein (1—n)Pu = QU —n), 


which is linear and can therefore be solved by the methods of 
Art. 21. 
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Example. Find the general solution of the equation 


dy 1 2,3 
Ge a> 


Divide by y’. 


dy 1 
—3 ee ee 
CG Fem Matas 
Let y? = u. oe 
pos 
ay dz dx 
du 2 ; 
avee 


t= — eS s* f we Sede + te) Pee aCe. 


Therefore — 22°y’? + cay? = 1 is the general solution of the 
equation. 

23. VARIABLES SEPARABLE 

Sometimes the equation Mdx 4+ Ndy = 0 can be brought to 
the form Xdx + Ydy = 0 where X is a function of 2 alone and 
Y is a function of y alone. In such a case the general solution 


is evidently 
Jf Xde + fYdy =, 


c being an arbitrary constant. 


Examp_e. Find the general solution of the equation 


avl — ydx + y~1 — 2dy = 0. 


Divide by V1 — 7? v1 — 2’. 


Therefore V1 — 2? 4+ v1 — y? =c is the general solution of the 


equation. 
3 
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The process of reducing the equation Mdzx 4+ Ndy = 0 to the 
form Xdx 4+ Ydy = 0 is called separation of the variables. 


24, EXACT DIFFERENTIAL EQUATIONS OF THE FIRST ORDER 
AND FIRST DEGREE 
Definition. The ordinary differential equation Mdx + Ndy=0 
where Mand WN are functions of x and y, is said to be exact when 
there is a function w(x, y) such that du = Mdz + Nady. 


Exampie. The equation 2«ydx + a*dy =0 is said to be 
exact because wu = 2’y is such that du = 2«ydx + a°dy. 

When there is a function w(x, y) such that du = Mdx + Ndy, 
then wu = c, where ¢ is an arbitrary constant, is the general solu- 
tion of the equation Mdx 4+ Ndy = 0. 

Condition that the equation Mdx 4+ Ndy = 0 be exact. If 
the equation Mdx + Ndy = 0 be exact, then, by definition, there 


is a function w(x, y) such that du = Mdx 4 Ndy. Now 
Ou 2 1OU 
du = ae dear a dy, 


from the definition of the differential of two independent variables. 


That the equation Mdz + Ndy = 0 be exact, it is therefore 
necessary that 
Cue 
Cy oe 
Conversely, the condition is sufficient. That is if 
oM _ oN 
Cy 0x . 
then Mdx + Ndy = 0 is an exact differential equation. 
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_aP 

* Ox = 
_@P eM aN 
“OyOx” Oy ~ Ox- 

man cor 3 / oP 

“On Oydr Ox\ dy 


Proof: Let ef. Mdz = P, 


Ore 
ir + Fy). 


oP oP 


where Q(y) is such that dQ(y) = F(y)dy. 


Therefore, if 
OM ON 


Oy Oz’ 


the left hand member of the equation Mdx + Ndy = 0 is an 
exact differential and therefore the aes is an exact differential 


equation. 
To find the general solution of the equation Mdx 4 Ndy = 0 


when the equation is exact. 
Let u(a, y) be a function whose differential is Mdx + Nady. 


ee TOU. 
Since an = M, 


u = { Mde + Fy). 
Since - = N, 
Bee, UN 5 Mda 4. — ioe. 


, dy) oO rs 


-, Fy) = i (3v- 5 if Med 
“.4u = falar + f (5 f atte) ay 


The general solution of the equation is u = ¢ where ¢ is an 
arbitrary constant. . 
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Exampie. Find the general solution of the equation 


(2° + 2ay + y)du + (y’ + @ + «)dy = 0. 
This is an exact differential equation. Therefore the general 
solution can be obtained by the above method. 


Since ae = M, 
Ox 


=f(@+ wy 4 y)de4 FY) =F + ey + ay + FQ). 


Ou 


O a 2 2 
Ei ov cee 


we pep vt ete 
, aF(y) . Pea ee 
co dy = y. FY) = 4 


y* 


a 
=_teytay4h 
y' 


Therefore ~ = ae xy + xy + =6 is the required general solution 
of the ea 


25. INTEGRATING FACTORS 
It sometimes happens that the differential equation 


Mdz + Nady = 0 


is not exact but becomes so when it is multiplied by some quan- 
tity. Thus, 


dy 
dg t Py = @ 


of Art. 21, is not exact but becomes so after multiplication by 


eS Paz, 7 
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Definition. A factor which changes a differential equation 
into an exact differential equation is called an integrating factor 
of the equation. 

Sometimes an integrating factor can be found by inspection. 


Exampe. Find the general solution of the equation 
(ae — y’)dx + 2xydy = 0. 


The equation is not exact as it stands but becomes so on multi- 
plication by 1/2’. 
Multiply by 1/2. 


: bie ie 2y 
a da + ~~ dy = 0. 
ede dee dy = 0. 
x x 


Aer y 
A . dé +a(L) =o. 


2 
Ss 47 <0, of ae Or, 

Therefore y? = — we* + ew is the general solution of the equa- 
tion. 

Rules have been devised for finding integrating factors in 
many cases where they cannot be found by inspection. For a 
discussion of them, the student is referred to Boole’s, Murray’s, 
or Johnson’s Differential Equations. 


26. EQUATIONS HOMOGENEOUS IN x AND y 


Definition. If Mand N of the equation Mdx + Ndy = 0 
are both of the same degree in z and y and are homogeneous, the 
equation is said to be homogeneous. 

To find the general solution of the equation Mdx + Ndy = 0 
when the equation is homogeneous. 


J ao ees 
dz = — N 
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i M ; 
Divide both numerator and denominator of — WV by « raised to 


the power indicated by the degree of M or N. 
Then every term in WM and J is constant or in the form of a 


coefficient multiplied by some power of = 


Then 

dy Toke 

ae (2): 
Let y = vx. 

dv 
. Ce + Vv = f(r). 
Therefore ; 
da dv 
Fe ORG 


an equation in which the variables are separated, and can there- 
fore usually be integrated without difficulty. 


Exampie. Find the general solution of the equation 
(2? + y’)dx — xydy = 0. 


dy v+y° 
dz vy 
Let y = vx. 
. gy AR nt C0 Le 1 
A Og as 2 eae a "Se ee 
dee 


.v = 2 log ex. 


Therefore y’ = 22” log cx is the general solution of the equation. 


27. EQUATIONS OF THE FORM ey aoa by +4 
daz ~ ae + by 1c, 
The general solution of an equation in the above form can be 
found as follows : 
Let =« +, and y=y' +%, where wz’ and y' are new 
variables, and a, and y, are constants. 
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Change the variables to z’ and y’. 


. dy ae + by +42 + by +4 1 
. dz! —_ / by’ b s ( ) 
4h a,% = 2 aie Xo at Yo =e C, 
Case I. If x, and y, can be determined such that 
a%,+5y+¢=90, and aa,+ by,+¢ = 0, 


then, on determining them such, equation (1) becomes 


dy! aya! + by 
Ce Soe by’ 


which is homogeneous and can be solved by the method of 
Art. 26. 


Case II. If, and y, cannot be determined such that - 
Cr Og a et: — 0, and dé, by + ¢ = 0, 
then, as was seen in algebra, 


1 


|.& 


Goes 
a b, ~ mM 


Ss 


2 
By substitution, the original equation becomes 
dy ae + by 4--¢, 


Oe maa + by) 4 ¢, 
Let a7 + by = v. 


d dy 
-& + b, “e = dx 
Therefore 
dv a+, 
7 coe Si 1inv + 6, 


an equation in which the variables are separable. 


ExamP.Le 1. Find the general solution of the equation 
d 
(2a 4 8y — 6) 7 = 6x — 2y — 7. 


dy 62,—2y—7 
de = I7'+. By — 6 
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Let c= 2,+2',¥=4,+ 4 where 
6x, —2y,-7=0, and 22, +4 dy, —6 = 0. 


dy Gx — 2’ 

. Ba G3 J as , . a 

tahoe, yal+y. OF gn eg ae 
,dv 6—2v 64 Sie 


i" dil = ooo . 
_ (24 80)dv_— dx 
6 dp 8 ae 

. — loge’ = log (3v? + 4v — 6)4. 


~ (3° 4 
say = Go + 40— 6h 


ee a ee at y—1 


Therefore 3y’ + 4ay — 62? — 12y + 14x = ¢ is the general solu- 
tion of the equation. 


Exampie 2. Find the general solution of the equation 


(82 —y +4) = 62 — 2y — 7. 
This comes under case II. 
Let 6a — 2y = v. 


-, v— 80 log (¥ + 38) = 2a +0, 
.°. 4x — 2y —30 log (6% — 2y + 38) = «, 


— ae a ae 
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Therefore 27 — y — 15 log (3x — y + 19) = ¢ is the general so- 
lution of the equation. 


EXERCISES 
Find the general solution of each of the thirty-six following 
equations. 
dy 1 ; dy 2 
i det gl=i-@. 2. dz + CObzY = cosec’ x. 
3. a =e — ¥, Ae Cleat a) _ +27y =a —2’ tans. 
dy x 1 
ee a 2). 
d 
6. 2(1 — at) 24 @ laze. = 
is wiles cos xy = $sin 2a. 
Sec wv ) ws - (2a? — 1)y = ax’. 
5 LL ? si TOG 1 yt xy = ary’ 
- 7, + siney = sine. - A — a!) 7 = ay = any". 
dy ae 2 dy ones 1 
jh de t C8 *Y =Y sin 22. 12. 3y it ¥ =*—!}. 
13. a — tanzy = y‘secx. 
14. y Vx? — Idx A, avy? — Idy 25 (Vy 
15. (& +1) cosvdz + & sin xdy = 0. 
16. ~2ay — y’ cosec a du + y tan xdy = 0. 
| d 
17. y(3 +9) 5 = 2(2y +8). 


i aay bei 7y" de 


+ (yf + 2y? — a 4 52*y — 2lay*)dy = 0. 
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19. (a + 4ay + y*)dae + (227 + 2ary + 47°) dy = 0. 
20. sin w cos ydx + cos x sin y dy = 0. 


21. (a? + log y)de + ; dy = 0. 


22. 2(@ — 2y)dy-- (a + 2y \de = 0. 

23. Saxydy — (2? + y’)dx = 0. 

24. (2 + 3ay — y’)dy — 3y'dx = 0. 

25. (a + 2axy)dy — (82° — 2ay + y°)dx = 0. 
26. dxydy — (40° + y’)dx = 0. 

20. (a —2ay dy + (@ — sey + 27 yd = 0, 
28, 8a°dy + (227 — 3y*) dz = 0. 


29. (Be + 2y — 7) 54 = Qn — By +6. 
30. (62 — by +4) Sh = ey 4. 
31. (5x — 2y +7) 5! = 0 By 42. 
32. (e— 8y 4.4) 2 = 5a — Fy, 

33. (a — 8y 4.4) = 20 — by +7. 
34. (5x — 2y +7) SY = 102 — 4y +6. 
35. (Qn — 2y 45) aay 48, 


36. (Gx — 4y 41) SE = 82 — 2 41. 


The following formulas, derived in almost any work in cal- 
culus, are inserted here for convenience of reference : 
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The subtangent and subnormal at a point (a, y) on a curve 
whose equation is expressed in rectangular codrdinates are 


dx d 
y= dy and y a respectively. The polar subtangent and polar sub- 


normal at a point (7, 8) on a curve are 7” = and = respectively. 


The angle between the radius vector to a point (r, 6) and the 
tangent line to the curve at the point is 


, 
dr” 
dé 

The equation of the tangent line to the curve y = f (x) at the 
point (2,, y,) on the curve is 


dy 


a a (4 — %,). 


“L=21 


The area enclosed between the curve y = f (x), the x-axis, and 
the ordinates whose abscissas are x, and #, respectively is 


[yt 


provided the curve does not cut the z-axis between x, and 4,. 
The length of the are of the curve y = f (2) between the points 
(&» Yo) and (x, y,) on the curve is 


ier 2) 


37. Determine the curve whose subtangent at a point on it is 
n-times the abscissa of the point. Find the particular curve that 
goes through the point (3, 4). Plot the curve (a), form = 1, 
(6), for n = 2. 

38. Determine the curve whose subtangent at a point on it 
is n-times the subnormal at the point. Find the particular curve 
that goes through the point (vn, 2). Plot the curve when 


nm = 4. 
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39. Determine the curve whose subtangent is constant and 
equal to a. Plot the curve, (a), when a = 1, (0), whena = 2. 

40. Determine the curve whose subnormal is constant and 
equal to a. Find the particular curve that goes through the 
point (1, 2). 

41. Determine the curve which is such that the length of the 
perpendicular from the foot of the ordinate of any point on the 
curve to the tangent line at that point is constant and equal to a. 
Determine the particular curve when e=a. At what angle 
does this curve cut the y-axis ? 

42. Determine the curve which is such that the area between 
the curve, the z-axis, and two ordinates, is equal to the are 
between the ordinates. 

48. Determine the curve which is such that the perpendicular 
from the origin upon any tangent line is equal to the abscissa of 
the point of contact. 

44, Determine the curve in which the angle between the radius 
vector and the tangent line is n-times the vectorial angle. Plot 
the curve when n = 3. 

45. Determine the curve in which the polar subnormal is pro- 
portional to the sine of the vectorial angle. 

46. Determine the curve in which the polar subtangent is pro- 
portional to the length of the radius vector. 

The equation for a circuit containing induction and resistance is 


( di : 
4 Lat Ri=e 


where ¢ is the electromotive force [E.M.F.] impressed upon the 
circuit, A the resistance offered by the circuit, Z the coefficient 
of induction, 7 the current, and ¢ the time during which the cir- 
cuit is in operation. In each of the four following exercises, 
determine the current in the circuit after a time ¢ supposing that 
the resistance and induction are constant. 

47. The E.M.F. is zero. Solve subject to the condition that 
@ = J when t = 0. 

48. The E.M.F. is constant and equal to E. 
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49. The E.M.F. is a simple sine function of the time, 
= Esin ot where EF is the maximum value of the impressed 
E.M.F., and o is the angular velocity, equivalent to 27n where 
m denotes the number of complete periods or alternations per 
second. 

50. The E.M.F. is the sum of the components each follow- 
ing the sine law, that is, e = E, sin of + E, sin (bot + 6). 

The equation for a circuit containing resistance and capacity is 


di a 1 de 
dawn .R at 
where ¢ is the E.M.F., F the resistance, C' the capacity, ¢ the 
current, and ¢ the time during which the current is in operation. 
In each of the two following exercises determine the current 
in the circuit after a time ¢, supposing that the resistance and 
capacity are constant. 
51. The E.M.F. is constant and equal to E. 
52. The E.M.F. is a simple sine function of the time, 
= Fsin of. 
The equation for a circuit containing resistance and capacity is 


where e is the E.M.F., P the resistance, C the capacity, q the. 
quantity of charge in the conductor, and ¢ the time during which 
the circuit is in operation. In each of the three following exer- 
cises determine the charge in the circuit after a time ¢, suppos- 
ing that the resistance and capacity are constant. 
53. The E.M.F. is zero. Solve subject to the condition that 
g =Q when t = 0. 
54. The E.M.F. is constant and equal to F. 
55. The E.M.F. is a simple sine function of the time, - 
= FE sin ot. 
ANSWERS 
z 


1. 4ay = 22? — x +. 2. ysin # = log tan 5 + ¢ 


38 


26. 


. cosec a + v2ay — y? — avers? 
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y=u—1+ce*. 4 Ye ek eee tee ee 


; x 
bs yl + x”) = log i+ (+2)! + C. 


; y= —5 log (aw akon, ie ge eee 


= ae + cx vt — 2 9. ; seh aeCly or 


Y) 
1 See 
Ae ae 


1 


2 
. f” = 2sin eg — a a elimsae abs 
ii 


Poet —— tee. 

y* = — 3 sin 2 cos’ & —sin’ & + ecos® x. 
ix? — 1 — seca + vy? — 1 — sec yoe. 
(%¥ +1) sine =e, 

iY 
a 
2y’ + by — 9 log (2y + 8) = 4a? +, 
a 


: Ree 
Z— Uy + Bey — Toy + B+ hy =e 


a4 
5 + 2ay+tay+tyt=e. 20. cosxcosy =e. 


e : oi 
: 3 tvlogy =«. 22, we? = e(@ +4)": 
3 
. (4y? — 2*)® = en?. ahs es 
(47 —#)° = ca 24. y (f=). 


2 
. ey + Bay — 827) = 2 {ou (G= 1)! 


2y + (8 + V21)e 


(Cg ee para oo 27. y= axlog <. 
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DS. Co = 


29. 
30. 


ol. 


32. 
33. 
34. 


by — (8 + v83)x] "3 
by — (3 — ion 


(y — 73) + 8@— Fs) — FD) —- (@-)’ =e. 
Oye =o) a C(4y — dy 3), 
{2(y — Bs) — (4 —V14)@ + 4)" (29? — Ay 


— 8xy + a’ + 28a 4 142) Ee 


ae gs) — (44 VI4)(2+4)}"™ 
(By — 5x 4 10)? = e(y —2 +1). 

15y — 30x 4+ ¢ = 3 log (5x — 15y + 17). 

dee — 2y + ¢ = 16 log (5x — 2y + 23). 


Lie Cm lor (ey 4.2). 
ZY? 16. — + log 22 — sy 4.1). 


L(ja+e 


2 


5 ali 


BOE 
36. 
37. y= ce; ¥ = FU. 
39. y* = ce’. 
40. y = 2ax + 
: Chie, swe 
Alle ya5(e+Se f 
42 an peli I haa) 
- Y= % (e =f 
44, 7” = esin n6. 
46. r = ce, 
R 
‘Bi = = 4+ce =" 
E 
49, t= 


) 


i ih 
38. Es tua ae Ss 


; yf = 2ax 4 4 — 2a. 


) y= 5 (e+ é *); zero. 
43. w+ y = ca. 


45. r=c— kcos 0. 


Bette 
AT. i=Te 4 


R 
(= sin wt — w cos ) +ce 2’, 


40 


50. 4 


51. 


52. 


53. 


5d. 


SHORT COURSE ON DIFFERENTIAL EQUATIONS 


E, aa 
+41= TR & sin wt — w cos ot) 
Laat o') 
Be 
+ R z e sin (bot +8) —bocos (bot +8) | 
L( i+ it 
R 
+ce & 
a 
4=—ce Fe 
i= ce Roy Eee (cos wt + ROw sin at). 
eel pe 
q= Qe Fe, 54. q= CE 4 ce #'” 


CE ~ aa 
q = 14 RC? (sin Obie RCo cos wt ) + COMRCS 


CHAPTER IV 


ORDINARY LINEAR DIFFERENTIAL EQUATIONS WITH 
CONSTANT COEFFICIENTS 


28. Definition. An ordinary linear differential equation is 
a differential equation in one dependent variable which is linear 
in the dependent variable and its derivatives. 

We saw in Art. 21 that the type of an ordinary linear differ- 
ential equation of the first order is 


where P and Q are functions of x, ‘and do not contain y or de- 


rivatives. 
In general, the type of an ordinary linear differential equa- 


tion is 


a" dy am 
that Paes coo bp Py a X 


and X, are functions of w, and do not 


wnere P,P, :-+, Pf, 
contain y or derivatives. 

In this chapter the only cases considered are those where 
P,, P,, +--+, P. are constants and real. Two forms of this equa- 
tion present themselves, namely, when the right hand member is 


zero, and when the right hand member is not zero. 
RIGHT HAND MEMBER ZERO 

29. We shall first prove a theorem used in the investigation 
of equations in this form. It is: 

Theorem. If y=y,, Y= %, *++, Y = Y,y are solutions of the 
equation 
dy dt ty d”y P 0 

a+ Pit Pat: Ps ni = 5) 
4 41 
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then 
Y= OY + Gs + + Oa GG ee ee 


arbitrary constants, is also a solution of the equation. 


Proof. Substitute y = ¢y, + ¢y, +-°:: + ¢,y, in the equa- 
tion. 


di” Om: Om 
ws a | et, wath, Tait + Pan] 


d™ am on 
a | get, ait Pit PES MS 5 Pa, | 


d"y, aa dy. 
+4) ge + igen ts Tait + Pate | = 0. 
Now each expression in brackets is zero, since 


Y= Up Y=Yy 2s Y=YMn> 


are solutions of the equation. Therefore 


Y = CY, + Ye + nay ae EE 
is a solution of the equation. 
Cor. Ify=cey,+ey,+°+:+¢,Y, isa solution of the equa- 


tion, then y = GY, Y = CYo °° +> Y = CnYn, are solutions of the 
equation. 

If y= 4, ¥ = Yy +++, Y= Y, are linearly independent solu- 
tions of the equation, then y ='¢y, + ¢y,+-°:-+6¢,y, is the 
general solution of the equation (see Art. 18). 


30. To find a solution of the equation 
dy ad” y d"™y 
aoe P, det + £; ice heat Fy =0; (1) 


in the form y = e””. 
Let y = e™ and substitute in the equation. If y= e™ is a 
solution of the equation, then 


em(m" + Pm 4 Pym'? 4 +.- 4 P) = 0. 
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Since e”* cannot be zero for any value of m, then must 
m" + Pm * 4+ Pm? 4.--4 P =—0. 


Therefore, if y = e™* is a solution of the equation, it is necessary 
that 
m” + Pm) + Pm? 4---4 P= 0. 


Conversely, if m has a value m, such that 
tHe Ae INP at ok Pe 


then y = e™*” is a solution of the equation. This is obvious be- 
cause on substitution of y = e™”, the equation reduces to 


ée™*(m,” ar Ej sf Fan a a ee aa PF.) =o 


Therefore the necessary and sufficient condition that equation (1) 
has a solution in the form y = e™ is that m be such that 


GPa bm” * 4. ee DP a O, 
Definition. The equation 
ee Hee HE ies oes tA =O 
is called the auxiliary equation of 


d”y (om ly a> *y 
ae PS + Poort: -+Py = 0. 
31. To find the general solution of the equation 


d"y Cay pty 
dx” fale seat P 2 dy aap Y 


weet Pof = 0. 


When the auxiliary equation has distinct roots. Denote the 
roots by m,, ™,, +++, m,. Then m linearly independent solu- 
tions of the equation are y = ¢™", y = e™*, +++, y = e”™*, and 
the general solution is y = ¢¢”"* + ee” + +++ + c,e™ (see 


Art. 29). 
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Exampe 1. Find the general solution of the equation 


2 
dy 30 aye 


dx? ° dx 


Let y =e. 
-,e™(m? —38m —4) =0. 2°. Cn —4) (m+ 1) = 0. 


Therefore y = ¢,e* + ¢,¢~ is the general solution of the equation. 
EXAMPLE 2. Find the general solution of the equation 


d’y 
dv? oo - 4y =. 
Let y =e 


e™(m? +m+4) = 0. 


(mt = ) (m4 1+ ws) =u; Seen 


1—V1y 14V15) ve : 
Therefore y =@e 2 *+4@e 2 * is the general solution of 


the equation. 
This solution may be written as 


vi5 iG 
y = ¢,e-** cos “ a+ ¢e7* sin = «% (see Art. 5). 


When the auxiliary equation has multiple roots. Suppose 
that the auxiliary equation m"4+Pm"7*4 Pm"? 4 ---4+P,=0 
has the roots m,, m 


m Mm 


Mag oth puiiee 

At first suppose that two roots are equal. Suppose for defi- 
niteness that m, = m,. Then a solution of the differential equa- 
tion is 


y = (G, ae oe, en" ae C,ense =e ¢,e"". 


Since ¢, + ¢, is equivalent to only one constant, this solution con- 
tains only n — 1 arbitrary constants and is not therefore the gen- 
eral solution of the equation. 

To find the general solution in this case : 
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Suppose that the differential equation is such that its auxiliary 
equation has the roots m,, m, + h, m,, +++, m,. The general 
solution of this equation is 


y= RU ais Cente oe Cer a aoc Cen 
= (SE ae Cer = Cane zs Se = Cer 


Expand e” by Maclaurin’s Theorem to » terms and the re- 
mainder. 


° iad Chay * Che \? 
fe es ee oes S fxs 
Malt hoy a + et be am, 
Oa ot 


Substitute in the above equation. 


-.¥ = (¢, + 6, )e"™ 
hee ee y* Os aly G 8 aa 
hz. myx hes Se ahieice ice 
+ eheere 1 4 4 rs eae ae 
| SG hil ta ect 
Since ¢, is arbitrary, h may be chosen such that ¢,h is any con- 
stant B for all values of h. Since ¢, is arbitrary, c, + ¢, may be 
chosen such that ¢,+ ¢,= A. Then 


hz hie ee ha)" 
y = Ae + Boom] 1 Sis iL = + ( ; “al 


Bs EE SRM oa (1) 
where A and B are arbitrary constants. 

Let fh approach zero. As fh approaches zero, the assumed 
auxiliary and differential equations approach identity with the 
given ones, and (1) approaches the general solution of the given 
differential equation. 


Now 
limit | 4% on (hx) (hx) ety 
[oa a n° a 


Therefore the general solution of the differential equation is 


Uy (A -+ ‘Br yen CCR A soir LG en”, 
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or, as we shall write it, 
Y= (¢ at Cu em ar Cer ap OOS Se cen 


In a similar manner it can be shown that if three roots of the 
auxiliary equation are equal, the general solution of the differen- 
tial equation is 


Y= (¢, + 6,0 + Ca em ae cee ae Ooo Se cen. 
and, in general, if 7 roots are equal, the general solution of the 
equation is 
ees GG a IC NES CG Ae eee 
Ifa pair of imaginary roots occur twice, the part of the general 
solution derived from these roots is 


v= (¢, ae CPN etal + (¢, 3 oe eke BL 
= (¢,+¢,7)e*(cos Bz +) sin Br) + (¢,+-¢,2)e* (cos Bx —j sin Bx) 
= €*[ (A + Bz) cos Ba + (A, + Bix) sin Bo]. 


Exampie 1. Find the general solution of the equation 


The auxiliary equation is m? 4+ 2m + 1=0, or (m4+1)?=0. 
The general solution is therefore y = e“(e, + ¢,). 


ExamMp_e 2. Find the general solution of the equation 


d‘y d*y dy dy 
i ae de ee 


The auxiliary equation is m* — 4m* + 8m? — 8m + 4 = 0, or 
Mi ne eg) ee 


The general solution is therefore 
y= &{ (A, + Bx) cose + (AL Be) sina}, 


32. As a physical application of the above principles, con- 
sider the following discussion (see Emptage, Electricity and 
Magnetism, page 180) : 
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In a galvanometer in which resistance is offered to the motion 
of the needle, the equation of motion of the needle for small 
oscillations may be written as 


dé dé 
dé + Qk di + w (6 — a) = 0, Cle] 


where @ is the angle through which the needle turns in the time 
t, kis a constant depending on the resistance offered to the mo- 
tion of the needle, w is a constant depending on the moments of 
the restoring forces on the needle, and a is the angle which the 
needle at rest makes with the line from which angles are meas- 
red. Let 6 —a = 6’, and substitute in (1). 
209) U 
i; oe £: 24 rh! =O, 

This is a linear differential equation of the second order with 
constant coefficients and right hand member zero. The auxiliary 
equation is m? + 2km + wo” = 0. ‘The roots of the auxiliary 
equation arem = —k+ 4? — ow. 

Casel. Ifk>o. 

In this case, 0— d= och’ 4. ¢.e(-#-V 4-09)! ig the 
general solution of (1). 

Case II. Ifk=o. 

In this case 6 —a = (¢,+ ,t)e&™ is the general solution of 
(1). : 

Case III. Ifk <a. 

In this case 0 a = &[c, cos Vw" — t+ ¢, sin Vo" — Ft] 
is the general solution of (1). 

In cases I and IL the motion is not oscillatory. The needle 
ean go through the position of equilibrium for one value of ¢, 
after which it reaches a position of maximum deflection and then 
continually approaches but never reaches the position of equi- 
librium. In case III there are oscillations in equal times, the 
periodic time being 
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EXERCISES 
Find the general solution of each of the following equations. 
Oy Uy ott Wo, = 
1. at ee St ae dit gq, — 2 =°9 
d°y d’ dy dy dy dy 
d*y d’y dy d*y d’y dy 
baa Saat ge le ee ee 
d®y dy dy 
d‘y d*y dy 
°: rin aos 5 Fa ar tae 
ANSWERS 
1 pet Aiea Sat 25% ee. 6 ena 
3. y= ¢e" + 66% + G. 4. y= e"*(4 + 6,2) + ce. 
Day = ee + xt ca’), 
Gli ¢.6° 40.6" *"( 6, COS ak eoaire). 
7. y= oe * + e*(e, cos % + ¢, sin 2). 
O. Yas CeO +L Cer" eL Cocos + e. Sin an 
3. y= (¢, + 62) cosa + (¢, + ¢@) sin &. 


RIGHT-HAND MEMBER NOT ZERO 
33. Symbolic form of equation. The equation, when the right 


hand member is not zero, is 


ac qu) pt *y 
a Po Pax =e + Py =X, (1) 


where P,, P,, --- P,, are constants, and X is a function of x but 
not of y. 
Let 
d 
Ste 5) 
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and the equation may be written, 
Dy + PD y+ PD y+ +++ 4+ Py =X. 


Suppose that y is treated as an algebraic factor of the left hand 
member of the equation. On this supposition, the equation 
becomes 


(Game DEP Dr. i. Py Se. 
Suppose also that D? 4 P,D”1+4 P,D"?* 4 .---4 P., factored 
as an algebraic expression in D, is 
CD ain) 9D) 22902); 
and that the equation is written 
(D0, )CD —m,) >=» CD — m,)y = X. (2) 


Equation (2) is not equivalent to equation (1) except in a 
symbolic sense. Let us see what conventions must be made in 
order that equation (2) be equivalent to equation (1). 

Let us make the convention that (D — m)u where u is a func- 
tion of x is equal to 


Also, let us agree that we shall begin at the right of the left 
hand member of (2) and work towards the left, evaluating 
according to the preceding convention at each step. Then 


dy 
(D <= m, )Y = dx et MY; 


(D —m,.,)(D —m,)y = (D —m,) (SE — ma) 


d’ d 
= a — (m,_, + ™,) + mY) 
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and finally, 
(D — m,)(D — m,) ++» D — my 


in n—1 

73 —(m,+m,4- * fm) see +774 (—1)?mym, -- my. 
Now —(M,+ M+ +++ +m) =P, +++; (—1)"mym, ++ +m i 
since the factors of D" + P.D** 4 P,D*? + .---+4 P,, treated 
as an algebraic expression in D are D—m,, D—m,, ---, D—*m,,. 
This expression is therefore the same as the left hand member of 
equation (1). Therefore, with these conventions, equations (2) 
is equivalent to equation (1). 


n= 


Exampie. With the above conventions, the equation 


d’y d S 
ae + Py Py=X 


may be written in the equivalent form (D—m,)(D—m,)y = X, 
where D—m, and D—m, are the factors of the expression 
D’ + PD + P, treated as an algebraic expression in D. For, 


(D—m,)y = oY _ mg, 


(D — m,)(D — m,)y = (D — m,) (2 = a 


d 
= 4 — (m+ m,) a+ M,mM,Y. 


Now — (m,+m,) = P,, and mym, = P,. Therefore the second 
form of the equation is eqinwalone to aN first. 


Definitions. When equations in the form (1) are expressed 
in the form (2), they are said to be expressed symbolically, or 
to be expressed by means of symbolic factors. 

When a symbolic factor D — m and a function wu are applied 
to each other so as to give (D — m)u or “ — mu, the function 

a 
u is said to be operated upon by D — m, or the factor D — m to 
be multiplied symbolically by w. 
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The factor D ~— m is called the symbolic Pperet or more 


briefly, the operator. 


34. Theorem. The order in which the symbolic factors in the 
equation of the last article are taken is immaterial. 

Consider in illustration the equation of the second order. 
Let the equation be taken in the form 


(D —m,) (D —m,)y = X. 
Then 


and 
(D — m)(D — m,)y = (D —m,) (5 — my) 
en ee 
Therefore (D — m,)(D —m,)y = X is equivalent to 
ts 
Also, (D — m,)(D — m,)y = X is equivalent to 


ay dy 
Fisk ok fan Py = X. (See Art. 33.) 


di 
nibs Py ad 


Therefore, in the case of the equation of the second order, the 
order in which the factors are taken is immaterial. 
The proof in the general case is left as an exercise to the 


student. 
35. First method of solution of the equation 


Let (D —m,)y =u. The equation then becomes 


du 
(D—m)u=X or dg = &. 
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The general solution of the equation 


du 
aes mu = X 
is (see Art. 21) 


us em i? e ™ Xda + ¢e™. 
.(D=-—m,)y=e"" hi Gar Xie ene 
y= em af emia if. ee XN dare 6c if dian 6 en 


— Ex Jf ear fem X(dx)? a eem ae Cen. 
This is the general solution of the given equation. 
Examp.e. Find the general solution of the equation 


a d 
y PUR Ph te 
ae oe + 2y = cos 2. 


Write the equation as (D — 1)(D — 2)y = cos2. 
Let (D — 2)y =u. The equation then becomes 


du 
(D—1)u=cose or =—w= cose. 


dx 
.uU= Bilas cos x dx + €,é 
= (sin x — cos 2) 4 Ge. 
.. (D— 2)y = (sin = cos 2) + Ge. 
*.y = $e" Hf @ “(sia — cos )deet Ce- cf é “dx +. ce" 
zp cos x — 8; sin # + ce" + ©,€. 
This is the general solution of the given equation. 
36. To solve the equation 
(D —m,)(D — m,)--- (D —m,)y = Xj 
we may proceed as follows : 


First, let (D —m,)-+-(D—m,)y = u. The equation then 
becomes (D — m,)u = X. 
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From this equation, w can be determined as in the case of the 
equation of the second order. Let 


(D — m,) ---(D—m,)y =». 
(D — m,)v = u. 


From this equation v can be determined in the same manner as 
uw was determined before. After n — 1 such steps there results 
(D — m,)y = 2 where z isa known function of z. 

The general solution of the equation (D—m,)y =z is the 
general solution of the original equation. 


Then 


37. The following theorems concerning the symbolic operator 
will now be established : 

Theorem I. A constant factor in a function may be written 
in front of the operator. 


Proof: Let au be a function containing a constant a as a fac- 
tor. Let D— m be the operator: Then 


(D — m)au = ee — mau, by definition 
y (3 mu) 
= dz oes 
= a(D —m)u. 


Theorem II. The result when the operator is applied to the 
sum of a number of functions is equal to the sum of the results 
found when the operator is applied to each of the functions 
separately. 

Proof: Let utv+w-t---+2be the sum of a number of 
functions. Let D — m be the operator. Then 


(D—m)(uU+vu-pw+t--> +2) 


ee nape eat ee) 
av 


du we mv ay nw sisee >= — Mze 
momen hiker et de OT ae 


=(D—m)u+ (D—m)v + (D—m)w+-++4+ (D—m)z. 


54. SHORT COURSE ON DIFFERENTIAL EQUATIONS 


38. The equation (D —m,) (D—m,)---(D—m,)y = X 
may be written in the form 


1 
aes (D= m, )\ CDi) TPG ITS te 


In the first form the symbolic operators 
D—m, D—m,, =, Dim, 


applied in succession give X. Moreover, by the theorem of 
Art. 34, the order in which the operators are applied is imma- 


terial. If the second form, therefore, is to be the same as the 
1 


(D—m,)(D—m,) ---(D—m, a 
must be such that, when operated upon by =~ 


D—m, D—m, ---, D—m, 


first, the symbolic expression 


in succession in any order, the result is X. 

1 
(D —m,)(D—m,) ---(D —m:; 
called the inverse symbolic operator, or, more briefly, the in- 
‘verse operator. 


39. Let 


Definition. The symbol is 


1 
heey aE A(T 


be a linear differential equation where the symbolic factors viewed 


as algebraic factors are distinct. Break 


uP 
(DS nes) ans) 
into partial fractions as if it were an algebraic expression in D. 


Then 


1 Brg i! 1 1 
. a — 1m, CD —m,) ~ i on Dene Dae a 
e 


: Mes a fil u XC = Oy. 


nm, —m, D— m, Be Nie 


Theorem. The result of operating on w+tv with 
CD, CD iii sical 


is 
A) 
4 YY 
ANS 


K 
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Proof: Operate on wu + v with (D — m,)(D — m,). 
(D — m,)(D — m,)(u + v) = (D—m,)(D— m,)u + 


(D — m,)(D — m,)v, by theorem II, Art. 37. 


Now 
(D—m,)u= A, and (D=m)v = — ae. X, 


— Mm, m, — m, 


Mm, 


by definition and theorem I, Art. 37. 
--. (D —m,)CD — m,) (u + v) 


iy ey, (D—m)(- J x) 


m, — ™, m — mM, 


1 


m,— Mm, 


[Cp ink — (D —m) Xx] =X. 


40. When the symbolic factors D — m, and D — m,, viewed 
as algebraic factors, are distinct, the result of operating on 


1 1 1 1 


y= = 
M, — Nt, — mM, m, —m, D'— m, 


with (D — m,)(D — m,) is X, by the preceding article, and the 
result of operating on 


vk 
1 (D — m,)(D — m,) < 


with the same factors is X, by definition. Therefore when the 
symbolic factors D — m, and D— m,, viewed as algebraic fac- 
tors, are distinct, the inverse operator of 


1 
Noms COE OCD = iK,5 


may be broken up into partial fractions the same as if it were an 
algebraic expression in D, and the result of operating with 
(D — m,)(D — m,) on the expression formed by multiplying 
each of the fractions symbolically by X, and taking the algebraic 
sum of the results, is X. 
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In general, when the symbolic factors 
D—™, Dm, <4 Dm, 
viewed as algebraic factors are distinct, the inverse operator of 


il 
’— Daj aD en 


can be broken up into partial fractions the same as if it were an 
algebraic expression in D, and the result of operating with 
(D — m,)(D —m,) ---(D—m,), on the expression formed by 
multiplying each fraction symbolically by X and taking the 
algebraic sum of the results is X. 

The proof of this theorem is left as an exercise to the student. 


41. Second method of solution of the equation 


i 
vote (D=—'m,)(D —m,) a 
1 
CD in) Da) 


if it were an algebraic expression in D. 


Break up into partial fractions the same as 


: 1 ete 1 1 
P<(D —2n, )\(D= m= tay n\n ie era 
Let 


1 1 i 1 
= D SOC Sil oS = Xe 
ih, Dean, nm, — m, D — m, 


Operate on wu with D — m,. 


du 
——— = hh 
re pi ae 
dx m, — ™m, 
1 
ss a CU Ce Kae oem. 
m, — mM, u 


Operate on v with D — m, 


dv 1 
A SO eS es SS 
. ere awad 
dx m, =m, 
1 
0 = = —— em ie da ee 
m,—™, = 
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1 my myx 
“.¥4¥= fag fe Adz — 


1 2 


Emme i emt Xx 


m, —m, 
an cm ae Ce, 
which is the general solution of the equation. 


Exampe. Find the general solution of the equation 


d’y d 


aA J 4 Oy = cose. 


dx 
Write the expression in the form 
a DAD) (D9) 
1 
=A (GOs 
it were an algebraic expression in D. 
J 1 ie 1 1 
Sree 2) D1 Do 


COS &. 


into partial fractions the same as if 


Break up (D 


Let 
1 1 
ao a and == con eV. 


Dt j Ne, 


Operate on wu with D — 1. 


.U=tcose —Fsin a + ce". 


Operate on v with D — 2. 


dv 
oi. 3 — 20 = — cose. 
dx 
‘v= 2e0s27+4 tain a + ce™. 


.°. ¥ = zy cos + — 58, sin x + ce” + ¢,e™, 
which is the general solution of the equation. 
This method does not apply when the symbolic factors viewed 
as algebraic factors are not distinct. 
5 
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42. It will be noticed in the example of the preceding article 
that the result is the same as that found by applying the method 
of Art. 85 to the same equation. This will be the case in any 
linear differential equation with constant coefficients to which 
both methods apply. 

The first method of solution will apply in all cases where the 
left hand member of the equation can be factored into linear fac- 
tors in D. The second method will also apply if the linear fac- 
tors in D are all distinct. If two or more factors are equal, and 
the inverse operator be broken up into partial fractions, the term 
or terms corresponding to these factors may be evaluated by the 
first method. 

Usually the second method is easier of application than the 
first. 


43. An examination of either method by which the general 
solution of a linear differential equation of the nth order with 
constant coefficients and second member not zero is derived shows 
immediately that the general solution consists of the sum of two 
parts, one containing terms not involving arbitrary constants, 
the other containing terms involving such constants. Moreover 
the arbitrary constants are involved so that when any one is zero, 
the term in which it appears vanishes. 


Definition, The part of the general solution of a linear dif- 
ferential equation with constant coefficients and second member 
not zero which contains the arbitrary constants is called the com- 
plementary function of the general solution of the equation. 


EXERCISES 


Find the general solution of each of the fourteen following 
equations. 


LA dy dy : d* d? d 
¢ JES rE / —s cae 4 y y y — 

- dx? T ee Ma eae oi 2 dx + 0 ae dz a. 
l4 dy dy | que eas : 

on dt ~ dy t 2Y =* 4, 5a 8 Kye 
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)_ d*y : d’y dy 
% Olay dit oqat’7=* 


6. oo 6S! gyi 
) 7. a en ae 8. Se eye 
0 9. C4 y = cosa. 10. Be ay we 
ff> 11. qetya sing. 12. oe - + Gf -y=cos-. 
/ 
13. ee 


In each of the six following exercises, find the equation of the 
elastic curve of the beam from the given differential equation, 
determining the constants of integration. Find also the deflec- 
tion of the beam. In these equations, Eis the modulus of elas- 
ticity, J is the moment of inertia of a cross section of the beam 
about a gravity axis in the section perpendicular to the applied 
forces, and / is the length of the beam. 


15. The beam rests on supports at its ends. It is supposed 
weightless with a weight P at its middle point. 
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16. The beam rests on supports at its ends. It is supposed to 
be of uniform cross section and of weight w per unit of length. 


17. The beam rests on supports at itsends. It is supposed to 
be of uniform cross section and of weight w per unit of length, 
and to have a weight P at its middle point. 


18. The beam is a cantilever fixed horizontally in the wall. 
It is supposed weightless with a weight P at its extremity. 


vp ea ge 


dx? 


19. The beam is a cantilever fixed horizontally in the wall. 
It is supposed to be of uniform cross section and of weight w per 
unit of length, and to have a weight P at its extremity. 


d’y 


1 


EI = Pie CE sete) 


WwW 2 2 
Stoll + 27). 
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20. The beam isvertical. It hasrounded 
ends. It is supposed weightless. It is de- 
flected a small amount a and a load P is 
applied at its upper end just sufficient to hold 
it in position. 


d’y 
EIT4 = Pay). 


The equation for a circuit containing re- 
sistance, induction and capacity in terms ot 
the current 7 is 


@i Rdi it Ide. 
ah Os. lade 


in terms of the quantity of charge’g is 
\ 


d’qg  Rdq q 1 
/ degli 10 - 1° 

where e denotes the E.M.F., R the resistance, Z the induction, 
C the capacity, and ¢ the time during which the circuit is in 
operation. In each of the three following exercises, determine 
the current and quantity of charge in the circuit after a time ¢, 
supposing that the resistance, induction and capacity are constant. 

21. The E.M.F. is equal to f(t). Solve when R’O x 41. 

22. The E.M.F. is constant and equal to FE. 

23. The E.M.F. is a simple sine function of the time, 
= Esinot. Solve when #°C x 4L. 


ANSWERS 
Bees ACO 06), 
er Sere Ce te 


I 
cole vd 
8, 


Eel Wey a = 
(24+-V 2 x (2—V 2 )x 
/ 3. Y= Sa + G,€ ; 
(8—V 13)x _ BV 1B)z 


ADY=edte ace 7 tee, 7 , 
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NT. 


18. 


19. 


20. 
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1_y? 5 me poe 
Y= pe? — Fer +4 +e “+ 46". 


ye al 
y= — a He AS Ce te 
y= (42 +. ¢, 4+ 6). Oe = © 42 eee), 
x sin & 


Me ee 


+ ¢, cos + ¢, sin %. 


fey ee A Ca ee cle. 


Xv COS & 


Y= + ¢, COs % +- ¢, sin &. 


2 


. ¥ = f(cos¢ —asins —xcosz) + ¢, sins 4+ ¢, cose 4 ce. 


ie 6 (Re Oar es +). 


y= — 2 — 242508 + oe” + esing + 6, cos &. 
dee W eee ‘ 1, Pe 

. AEly = pe as Deflection = Er 
i ka TE : 5 wilt 

. SEL = VU gt Deflection = 384 ET 


WE SEP 2) 90 hehe Pye wl* 
Gor 1 See Ge eerie 
LORE 5 wilt 


OETy = — 


Detection mr ean? 
2EIy = — Pla’ + ao Deflection = ee 
2ETy one ee eee ; 
eo Deflection = ae Salk spy! ts 
v= oe vers z. 
Deflection = a vers ee and.°, P= ake 
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C : is a ee 
2is= ee Ff eFpctyat =e 2 {Apa 


where 


eee ands ie ihe 


RF? 0? — 410 PR Cua CoA Cs 


The value for q differs from that for 7 only in having f(t) instead 
of f’(¢). 
_ RC-NR2O2-ALC , _ RCN R?2C2-ALC , 

22. 1= Ge aoe! + ¢,¢ eal z 
when f?C > 41. 

—3t VALCO — RC? —gt. VIL — FO 
Be 2L TAS 
1= ce cos ILC acs sin oLO ; 
when #°C < 4L. 


a: Rk t _£ t 
ta 6.6 Pete ahs 
when R’?C = 4L. 
RC-JR202-410 _ RC+ JR202-4L0, 
2LC0 


q=CE+ce UC +e¢ 


when R?C > 4L. 


-%t V4LC— FC? 
2L vin 
q= CE + ce cos oLC t 
+ ¢ pth gee 
: 2LC ; 
when f?C < 42. 
_-, maeie 
q= CEL ce + efe ~, 
when A’C = 4L. ; 
Bol ig fe 
2 
23. i= ioe ; Sin wt + e , COS wt 


Beets) 24(4 1) 


t t 
+oe T4+e¢ %, 
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where 7, and T, have the values given in exercise 21. 


a(t 


ar (1 i) 


— HwoR 
q= 1 
Pw? + (z — Le! 


, COS wt + ; Sin wt 


t 
-;% = 
+¢e 14 6,6 


t 
Tp. 


CHAPTER V 


HOMOGENEOUS LINEAR DIFFERENTIAL EQUATIONS. EXACT 
LINEAR DIFFERENTIAL EQUATIONS 


HOMOGENEOUS LINEAR DIFFERENTIAL EQUATIONS 


44. Definition. A homogeneous linear differential equation is 
an equation of the form 


n, n—1 


n'y a dy d 
a et PO et to + Peat get Py = X 


where P,, P °**> Pn» P, are constants, and X is a function of 
x but not of y. 

This equation can be (oie reece into an ordinary linear dif- 
ferential equation with constant coefficients by changing the inde- 
pendent variable from x to z, the equation of transformation be- 
ing z=e. The equation that results from the transformation 
may be solved by the methods of the last chapter. If a solu- 
tion is y = f(z), the corresponding solution of the original equa- 
tion is y = f(log x). 

45. The transformation and general solution of a homogene- 
ous linear differential equation in the general case will not be 
considered here. We shall merely consider them in a particular 
example. 

ExamPpLe. Find the general solution of the equation 


Dy d’y dy 

3 ee eh aloe 
dx? t ie dx? : da Ne: 
Let « = & z= log a, 


dy dydz 1 dy 

dx dzdx x dz’ 

dy d € Z) eden Cisdey dey ot (3 | 
~ dx he rc : 


2 dz) 


AD Ema Pon Tee meee 
See Gs. 
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Sul aided i ee LY I Te 894 +24). 
dF me da Vat Nz? dz x \ dé dz dz 


dy _ dy 
da~ dz 
vy dy dy 
da?” dz? dz 
d’y 


The general solution of this equation can be found by the 
methods of the last chapter. It is 


y= — fz aF t ar ce” == tec aF ce". 


The general solution of the original equation is therefore 
Gq C, 2 
y= —thget hp 34 Sy ge 


EXACT LINEAR DIFFERENTIAL EQUATIONS 
46. Definition. A linear differential equation 


d” Oe d* d’y 
Je at P Je sa op Pa Pas rath, n—-1 mg x 
da da dx dx dx 


is said to be exact, when, if the left hand member be represented 
by V, the expression Vdz is the exact differential of some func- 
tion U which does not contain an integral of y. 

The expression U is evidently an expression actually contain- 
ing a derivative of order n — 1. 


47. To find the necessary and sufficient condition that the 
equation of the preceding article be exact, and a method of solu- 
tion of such an equation. 
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Multiply each term by dz and take the integral of each term. 


at a eer dep f iA 
pts aaies {Peg de + { Pyde = [ Xie +. 


; ie Pydy = if Pydy identically. 


And, by integration by parts, 
(ee = dz = = = ( Pi_dae or Peay 
d’y Y 
Page 3A = a [P ni 1x —P anid te Tea 
d* Yy vty | YA / / ” 
fPuagd g de = = as 1 nse +— ip n—3 ink Z n—3d =F Nes) ? 


where the primes denote differentiation with respect to x. 


Pape = pee es ae Pe = Pi ee yd 


Se ee eg A ke OY 
fer eet ee ie 
tyson eee 
ae 

Write the expression in brackets as Q,, 9,1, °+*, Qo respec- 


tively. 
q7 
(de pom f Qyde+ Qt Waite + Qype C1) 


Now in order that the equation be integrable there must be no 
term in the right hand member of (1) containing an integral of 
_y. The necessary and sufficient condition for this is that Q, =0. 
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Therefore the necessary and sufficient condition that the equation 
be exact is that Q, = 0 or 
Ey in) A cea, si bi al ee 5 Si 0. (2) 
When this condition is satisfied the equation reduces to 
ary dy 
Qt to HOt Qa = fAdepe (8) 


Tf the coefficients in (8) satisfy a relation in Q similar to (2) 
in P, equation (3) is exact and the above process may be re- 
peated. 


ExampLe. Solve the equation 
ay ay dy | 
2 “eee pet — 
ei a) a8 + (1 — 52*) da Cie + 2y = 6x. 
Here Po: 2).— Der= 2,2 af Oe 10 
ae eaten gle 3 EU Sle Tee) ic er | 


The necessary and sufficient condition that the equation be exact 
is therefore satisfied. 


Oa — 274+ 10t Oe] 2, Oli be 1 oe = 


eee Os 
The equation therefore reduces to 
gee 
In this equation, P, = 2%, — P’,_, = 4a, P”,_, = — 62. 


iP oP eae 


n—1 
The necessary and sufficient condition that this equation be exact 
is therefore satisfied. 


0,4 = — 24327 —Llee’—l, °Q, sax —@), 
The equation therefore reduces to 


o(1 — at) 24 @— Ly =o + 00 + ey 
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This equation is not exact. It is however an ordinary linear dif- 
ferential equation of the first order, and can therefore be solved 
by the method of Art. 21. The general solution is 


2 


g= — 40 lor (1 — 2) + fax log —~— a 
nee + 4elogs +2) + ¢f, 


which is therefore the general solution of the original equation. 


EXERCISES 


Find the general solution of each of the following equations. 


i BTS + Tega + Be! = (log 2) 


2. 2 ig aya 

3 ool eT 4 6x9! — by =a? 

4 Nem ee 4x “Y 4 Ay = log a 

5 oot 4 tt ey 0 

6 po =e if oot oD yy loge 


2 
8. Sat eet y-™ 


Sh ee 5 do 4 Dy =o. 


d’y 


dy 2 
10. at Sa + cosec? ry = COS 2. 


ine (@t—2)§ Ge 2 +y=0. 
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d’y dy é 
22 eae — 
12. (+ 32") Aah 2d + 627) pie 6y = sin x. 
a & 
13, (4 o—8e4 1) 4 Or Oe 8) 8 


d 
+ (18x 46) 5/4 by =a". 


dy  , dy | dy 

ope el - niche ee 

14. x dt 5x dat + as =— 
ANSWERS 


1. y = t(log x)* — $(log x)’ + $$ loge +6, 4 = aE e 
¢, 
2. y= — Bt + Oe + oe 
3-4 = 42 lor a — feat ee 4 ee 4+ 62°. 
4, y= floge4 dt ee + qa" 4 4. 
Gn eos c 
5. y= t+; loge + qe. 6. y= Re pote, 


7. y = loge + 2 4+ ¢2 logs + on. 


x .% ¢ 


a: v 2 
9. Ue ge ale or setae 


10. y= “sin x + ¢, sin x log (cosec x — cot x) 4 ¢, sin 2. 
Sena cs & 
ll. y= Ee lor(@:— 1) 5 re 


12, (@+ 32’)y = —sing+e¢r+46, 
13. (a2? 4 a? — 84 4. 1)y = gh +6, + 60 4 62%. 


1 
14. y= 5, (log x)? +o, + log a + +. 


CHAPTER VI 
CERTAIN PARTICULAR FORMS OF EQUATIONS 


48. An equation in the form S —<, = f(2). 


An equation in this form is exact and can therefore be inte- 
grated by the methods of the preceding chapter. It can also be 
integrated by direct integration. 

The first integration gives 


gd? 
Tent = J Made + 4 


where a, is an arbitrary constant. 
The second gives 


F=f [ide + ae + 04 


where a, is an arbitrary constant. 
After n integrations there results 


g= fff SIO Gey pee opt + +e, 
where ¢,, ¢, +-+,¢, are es constants. 


49. An equation in the form & —~ = f(y). 


An equation in this form can in See be integrated only 
when n = 1 and n = 2. 
When n = 1 the equation is 
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To integrate, separate the variables. 


When n = 2 the equation is in the form 


= f(y). 


To integrate, multiply by 2 


dy d’y @ f dyn? 
dx dx? dx ; 


“ae (GE) = 27 SE. 
(G1) = 2fFaddy +4, 


Suppose that'2 { f(y)dy = W(y). 


d a 
peed 


dy 
Siemens ee, 
dy 
et+e¢=2+ | ——— 
ee a Gs 


50. An equation that does not contain x directly. 
Such an equation is of the form 


dy d"y 
F(y¢ ’ + aS) = 
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Let 
dy 
dpe 
Then 
ay dp dpdy dp | 
da dx = dydz = Pay? 


dy d(/dy\ d/ dp\  d/ dp\dy 
da? = da \ dat) = de a) =a P ay} dx 


and so on. 

The equation then becomes a differential equation in p and y 
of order n — 1. Suppose that it can be solved and that the solu- 
tion isp = f(y). Then a solution of the original equation is 


ion ee 
has GY 
51. An equation which does not contain y directly. 
Such an equation is of the form 


F( Lo Tn) =O 


fp et 
> dx’ > dx” 


Let 


The equation then becomes a differential equation in p and x 
of order n — 1. If the equation can be solved for p and the 
solution is p = f(x), a solution of the original equation is 


y+eo= f f(a)de. 
52. An equation of the first order solvable for y. 


In such a case, when solved for y, the equation becomes 


y = F(, p). (1) 
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Differentiate with respect to x. 
d 
ey Oe (x P; 2). 


This equation does not contain y explicitly. It is an equation of 
the first order in p and x. If it can be integrated as an equation 
in p and 2, there results on integration an equation between 2, p 
and an arbitrary constant. From the resulting equation and 
(1), if p can be eliminated, there results an equation between 
x, y and an arbitrary constant, which will be the general solution 
of the equation. 


53. An equation of the first order solvable for x. 


In such a case, when solved for x, the equation becomes 
y= Fy, p). Differentiate with respect to y. 


1 a 
a Y, (oa he 
P (2 dy 


The method of procedure from this point is similar to that in the 
preceding article. 
EXERCISES 


Find the general solution of each of the twelve following 
equations. 


ile ae 2. CY = cos 
3 oF — loge, 4, . — ay. 
5 HL eae, 
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d’y dy dy dy 

9. © Tat 8a = 8a. 10. ¢@= 7 +o 
dy \* dy (dy\’ 
i. say + (3). 12, ya aGt (3). 


13. Find the curve whose curvature is constant and equal 
to x. 

14. If a sphere of radivs A, is surrounded by a concentric 
shell of radii #, and F,, the potential function, V, at a point 
either in the space between the conductors or outside the outer, 
satisfies the equation 


PV 2dV 0 
dea Oro 
where r is the distance of the point from the center of the sphere. 
Solve the equation given that V, is the potential on the sphere 
and JV, on the spherical shell. 
15. If a circular cylinder of radius F, is surrounded by a cir- 
cular cylindrical shell of radii A, and R,, both of very great 
length, the potential function, V, in the space between the con- 


ductors, is such that 


a View oy 


dr’ r dr 


0, 


where 7 is the distance from the point to the axis of the cylinder. 
Solve the equation given that V, is the potential on the cylin- 
der and JV, on the spherical shell. 


ANSWERS 
see te Oe bret EG 


antl 


eT 


Ce a 


2. ¥ = —coS7+¢,%-+ ¢,. 
3. y= — log « — +(log x)? + ¢# + ©, 
4. y=¢,sin(ar 4 ¢,). 
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14. 


15. 
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+B+ 6, -¥ Y ey — 2 += “log (vey + Ve — 2). 


C120 
Arm Cle. 


+2+0= ; log (cy + vey? — 1) 


go oo abe ae) nee 
1 oY 

y= atte, 10-9 2 Ae eee, 
dete = eve —y — log (1 = ve — 9). 
y= ce —C. 13. A circle of radius = 
Ve RR, V,—V, , &,V,—R,V, 

Spee Cee Re 

Vm Pog e 4 EAT te, 


log R. log —} R, 
2 


CHAPTER VII 


ORDINARY DIFFERENTIAL EQUATIONS IN TWO DEPENDENT 
VARIABLES 


54. So far, the differential equations considered consisted of 
two variables, one independent and one dependent. We shall 
now consider equations in three variables. These may be divided 
into two classes : those in which there is only one independent 
variable, and those in which there is only one dependent vari- 
able. The first comes under the class called ordinary or total 
differential equations: the second, partial differential equations. 
This chapter is taken up with a discussion of a few forms of ordi- 
nary differential equations. The next chapter is devoted to 
partial differentia] equations. 


55. If f(x, y) isa single valued and continuous function of 
the two independent variables x and y, given by the equation 


Oz Oz ; 90 
z= f(a, y), and a, and ei continuous, then, by definition, 


Oz Oz 
dz = ay = By 


or 


of(, y) of ») 
dz = an da Bu dy. (1) 


If f(x, y, 2) is a single valued and continuous function of the 
three independent variables w, y and z, given by the equation 


Ou Ou Ou : 
= fi 
u= f(x, y, 2), and ae Oy and Dg are continuous, then, by de 


nition, 
Of (a, Ys #) Of (&, Y 2) Of (x, Y, 2) 


77 
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56. Equation (1) of the preceding article has, as a special 
case when z = 0, the equation 


LOD ay 4 ODay — o, 
oy 


That is, the equation is true for the equation f(z, y) = 0 where 
x and y are independent variables. If y in f(z, y)=0 is a 
single valued and continuous function of x, the equation holds 
true for all values of 2 for which y is a single valued and con- 
tinuous function, for in this case y is merely restricted to values 
which it could assume, as well as others, in the more general 
case where it is independent. 

This can be seen more clearly perhaps by a consideration of 
the geometrical representations of the equations. 

The equation z = f(x, y) when x and y are independent vari- 
ables represents a surface. Ifz = 0, the surface is the xy-plane, 
and the equation 

we ee vs oe y) tear 


holds true for every point in the plane. If y is a single valued 
and continuous function of a, the equation f(x, y) = 0 repre- 
sents a curve in the xy-plane in which the equation expressed in 
the form y = $(#) gives a single valued and continuous function 
of x, and since 


f% Y4 , KY) 
rae de + by dy = 0 


holds true for all sets of values of x and y in the plane, it holds 
true for all sets of values which together determine a point on 
the curve in the plane. 


57. Equation (2) of Art. 55 has as a special case when u=0, 
the equation 


— Be 2) ap 1 SG H ® ues _ 2) 
dzi4- By dy + dz = 0. 
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By reasoning similar to that employed in the preceding article 
in the case of two dependent variables, it may be seen that 
this equation holds true when z is a single valued and continuous 
function of and y. 


58. An integral relation in x, y and z, equated to an arbitrary 
constant ¢, say $(2, y, z) = ¢, can always be expressed in the 
form 


Pax + Qdy + Rdz = 0, 


where P, Q and F are functions of z, y and z, and do not con- 
tain the arbitrary constant c. 

For, the result of taking the differential of each member of 
the equation $(2, y, 2) = cis, by the preceding article, 


O(%, Y, 2) OP(%, Y, 2) O(a, Y, 2) 


and this equation is in the specified form. 


Examp.e. The result of taking the differential of each mem- 
ber of the equation 2’y — az’ + y’z = ¢ where ¢ is arbitrary, is 


(Qey — 2) dx + (27 + Qyz)dy + (4? — 2xz)dz = 0. 
_ This equation is in the form Pdx + Qdy + Rdz = 0. 


| The resulting equation Pdx + Qdy + Rdz = 0 is such that P, 
- Qand F£ are proportional to 


| Op O¢ Ob 


=, = an 
Ou’ Oy Oz 

respectively. 

| Conversely, however, an equation of the form 


Pda + Qdy + Rdz = 0 


where P, Q and FR are functions of w, y and z, does not neces- 
| sarily give rise to a solution of the form $(a, y, 2) =c¢. This 
| can be seen immediately because an equation of the form 


Pde + Qdy + Rdz = 0 
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which gives rise to a relation ¢(a, y, 2) = ¢ must be such that 
P, Q and £ are proportional to 


ap op |, oO 


Ox’ Oy oe es 
respectively, and these relations cannot hold for all values of P, 
Q and R. 
59. To determine when an equation of the form 
Pdx + Qdy + Rdz = 0 


has a solution of the form $(2, y, 2) = ¢. 
If it be assumed that Pdx + Qdy + Rdz = 0 has a solution 
p(x, y, 2) =¢, then P, Q and F must be proportional to 


respectively, or 


where p» is a certain unknown function. From the first two of 
these equations there results 


or 


Oz Cy’ : Ce 
Similarly, by using the first and third equations we get 
ORM CE Op. > Om 
(5-5) PERE (2) 


and by using the second and third, 


(5-5 ) = Re of. (3) 


Oz «Oy 
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Multiply equations (1), (2) and (3) by R, Q and P respec- 
tively, and add. 
: 0Q OF Che cl oP oq 
+ P(E 5) + (5 — + ae )= 2 
Therefore, if the equation Pdz + Qdy + Rdz = 0 has a solution 
$(4, y, 2) = ¢, equation (4) must be satisfied. 
Conversely, if equation (4) is satisfied, the equation 


Pdz + Qdy + Rdz = 0 


has a solution (2, y, 2) =¢. The proof of this theorem is some- 
what long and will not be given in this book.* The theorem 
however will be assumed in the subsequent work. 


Definition. Equation (4) is called the condition of integra- 
bility of the equation Pdx + Qdy + Rdz = 0. 


60. To solve the equation Pdx + Qdy + Rdz = 0 when the 
condition of integrability is satisfied. 

Suppose at first that z is constant so that the equation becomes 
Pdx + Qdy = 0. Solve this equation. Suppose that the solu- 
tion is f(x, y, 2) =a constant. Let u=f(a, y, 2). Find a 
quantity » such that 


Ou 
pP= Ox 
Ou 


Multiply the equation Pdz + Qdy + Rdz = 0 by p. 
.°. p Pda + Ody + Rdz) = 0. 


This equation may be written in the form du + Sdz = 0 where 
u and S are in general functions of 2, y and z. In the equation 
du + Sdz, change the variables from a, y and z to x, wu and z by 
means of the relation wu = f(%, y, 2). The equation then be- 


* For a proof of this theorem and also that S’ of Art. 60 does not con- 
tain x, the student is referred to Forsyth, A Treatise on Differential Equa- 
tions, Art. 152. 
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comes du + S’dz = 0. It can be shown that S’ does not contain 
x. Assuming that it does not, the equation du + S’dz can be 
integrated as an equation in wu and z. The general solution of . 
the equation is the general solution of the original equation. 


ExaAmpLe. Solve the equation 
YZ ae . UZ 
oppo wy 


Suppose that z is constant. The equation then becomes 


{dy — tan” dz = 0. 
x 


We ie 
4y oo 


The solution of this equation is 


=—,dy=0 or yde—xdy=0. 


x 
7 = a constant. 


Let u = —. 
y 
ou 1 
Ox = y 
yee 
S y 2 2 
Wey 
as 
: ee : ge aly 
Multiply the original equation by —* 
yz 
2 2 
can ee =s * dy = TTY tant dy = 0. 
y y i x 
Now 
iL z 
du = —dx — = dy. 
y got 
2 2 
Sites TTY tan g, = 0. 
z z 
Substitute 
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in this equation, y being derived from the equation 


© 
Lm) er, 
y 
ne i UE ds = 0. 
Z u 
Separate the variables. 
du dz 


1 
2 1 
Cr +1) tan™— 


dv Zz 
a = ()) 

v ap Zz 

Vz = C. 
ztan ? y IC! 
oF 
Therefore 
z tan? J = C 
x 


is the general solution of the original equation, 
61. Suppose that in the equation 


Pdx + Qdy + Rdz = 0 
the condition of integrability is not satisfied. Then there is no 
relation ¢(z, y, 2) = ¢ which satisfies the equation. In such a 
case a relation 
W(x, y, 2) = 90 

is assumed arbitrarily and a relation $(2, y, z) =e is sought 
which, together with y(2, y, z) = 0, will satisfy the equation. 
By differentiation of y(x, y, z) = 0 there results 


oY de +5 = 0. 


From this equation and (1) suppose that z and dz be eliminated. 
Then there will result an equation of the form P’dx + Q’dy = 0 
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where P’ and Q’ are functions of « and y the values of which 
depend upon ¢(a, y, 2). Suppose that a solution of this equation 
containing an arbitrary constant is found and is ¢$(2, y, z) =. 
Then this solution and y(2, y, z) = 0 together give a solution 
of the equation. 

As an illustration consider the following example : 

The equation 


LG + Kine, 


considered in exercises 47 to 50 inclusive, Chapter ITI, for special 
cases of ¢, does not satisfy the condition of integrability if e, 7 
and t are variables independent of each other. For, the equa- 
tion may be written as 


Ldi + (Ri —e)dt + 0-de = 0. 


By application of the condition of integrability there results 


(079 Sl eee 


yl Cl, (OCRN =e me 
+0{S— oi } =o, 


or 
ot zai); 


Since LZ is not zero, the equation does not satisfy the condition 

of integrability. Assume e = f(t), however, and the equation. 
becomes an ordinary linear differential equation of the first order. 

The solution is 


at fe eae 
i=> f et fat + ce x 


From this solution the results of exercises 47 to 50 inclusive, 
Chapter IIT, may be found by substitution. 


62. The cases considered thus far consisted of one equation in 
two dependent variables. Another important class of equations 
is the case of two total differential equations in two dependent 
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variables where each equation is of the first degree with constant 
coefficients. The method of solution of this class of equations is 
as follows : 

By differentiation and elimination, obtain one equation in one 
unknown. ‘This equation may be solved by methods previously 
discussed. The solution found must be a solution of the original 
equations. Another solution is found by substituting the one 
just found in the equations. The complete solution consists of 
two linearly independent relations between the variables. 


EXAMPLE. Solve the equations 


dz 
de 7 J +2 = 29: (1) 
dy 


Differentiate (2) with respect to x. 
ty ody dz 


73-65 +57 = 0. (3) 
Multiply (1) by — 5, and add to (2) and (3). 
dy . x dy 


gg OG + ty = 9. 


This is a linear differential equation of the second order with con- 
stant coefficients and right hand member zero. It can therefore 
be solved by the methods of Art. 31. 


Sp tee otto (4) 
Substitute this value of y in (2) and solve for z. 


ea Oe 4 & of, (5) 


Equations (4) and (5) together constitute a set of solutions of 
the given equations. 
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EXERCISES 


In each of the seven following equations, show that the condi- 
tion of integrability is satisfied. Solve the equation. 


1. Y+e2)de+ @+a)dy+ (@+y)dz = 0. 
2 Ory 2a i oye se dae (ae ey ae = aie 

+ 2y'z + Qy2? + 1)dy+ Cy +¥+y'2+4+1)dz=0 
(Qay + 2*)dau + (a? 4 2yz)dy + (y’ + 2az)dz = 0. 
(a + 2)ydx + (a + 2)ady — xydz = 0. 
(y + b)(24+h)dx+ (ua) (2+h)dy+ (+a) (y+6)dz=0. 
(ye + Qn) dau 4 (a2 + Qy)dy + (ay + 22)dz = 0. 
(Qaya + ya y2?)dx + (a%2 + 2xyz + wz’) dy 

+ (ay + xy’ + 2xyz)dz = 0. 


Solve the following sets of equations. 


ee a tee 


dz 7 ty 


dz ody dy 


d dz 
10. + 3y +2 = 0, o, + 38y + 52 = 0. 


dy dz - 

11. 7 + 8y + 22 = 0, dat *¥— 4 = 
dy dz 

12. de 7 oY — 2% = 9, det ¥ — 2% = 9. 


dy daz 


dz da d 
tq tby—saete, = 


15. 


16. 


17. 


18. 


19. 


IDA 
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dz d 
det Yt 8% =% vt fy + eae. 


dz d 
de 7 OY + = % tty eae. 


dz d 
de + OY — 2% =% at WteHe. 


dz e dy 2x 
de t ¥ — = % = ae 
dy 32 dz 4a 
de 7 ey 2 = det OY — 22 =e. 
ANSWERS 


ry + Y24+ 2% = C. 


. wy + y'2 + log(a+y+2)=¢6 

. Sy tye+ ve = ¢. 4. cy=c(a+2). 
Pitt OY Oya ket, 6. x72 bia 1 4? oe © 
. eyelet y +2) =o. 


7 7¢ Té 
—8a. — 32 1 23a 2 7-32 Sip ee. 
- Y= Ce + ¢,7é ; a ee ahaa: 2 Heng 4 
8¢ 
2 —6x —x 2 5-62 
y = C,é + C,€ ’ | c,é — a é 5 
—6z —2a — tae he 
-y=oe”*+ee"%, z= dee" — ce 
1416 , 1+ 6, 
y=ce 7 +46, ; 


1+ 65 


res x ¢, 
2=— 3 (7 + Vo5)e pe (7 NOB je. 2 


14 65 


Ee ee le Bei eg aes 
y= + ¢, ) 


2 2 


suk = ON ead OR ail 
e=( 7h ~$)¢ cos ye — (f+ if Je OE eae 
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13. y = ¢e* cos Vox + ¢,¢7 sin y5a, 

Z= —C, v5e™ sin V5a + ¢, V5e cos V5. 
oy = — 2% — tte ee - ee, 

@=n — 1 — 2x — $e* — Exe?’ 4 Bee + ,e~. 

= C2 et 

15. y= pe" Fe +¢6° +66", z= = —re"+he— Se "456 ™. 
16. y= pe ale he? al: C, e7@-¥ 3) 3)x ES C, Pecans ae 

a ae" +34 le ae C, V3 Can V 3)z ie. 3 eat 3a 
LT. y= — rhe tae +e 4 Oe 


t= — tes + ote (8 4 VID 
ia ¢,(3 = 14) “Ot e, 


18. y = 4e* — fe" + ee” cos Vor + ¢,e* sin dz. 
2= — ter 4 (3 he | cos Va 
4 i= 
“be ae Nde, e* sin Vda. 
3 3 
19. y = pre* + te 4 ce cos V10x + c,e* sin V10z. 


a= = aye felt 08 ot 005 10: sn: 


sin ¥102. 


CHAPTER VIII 


PARTIAL DIFFERENTIAL EQUATIONS. 


63. So far we have considered differential equations in which 
there is only one independent variable. We shall now consider 
equations involving two independent variables. Such equations 
belong to the class called partial differential equations, 

In this book, the independent variables will be denoted by « 
and y, and the dependent variable by z. ‘The partial derivative 
of z with respect to x and with respect to y will be denoted by p 
and q, respectively. 


Definition. A linear partial differential equation of the first 
order is an equation of the form 


ie Pt Q4q a7 R, 
where J, Q and F are functions of 2, y and z, and do not con- 
tain p or q. 

64. If there are two equations containing x, y and z, p and q, 
which can be solved for p and q, the result may be substituted in 
dz = pdx + qdy 
thus giving an ordinary differential equation. Usually, however, 

there is only a single differential equation given. 

65. Derivation of a partial differential equation. 

(a) By the elimination of constants. Let $(a, y, 2, ¢, ¢) =90 
be a relation between x, y, z and two arbitrary constants c, and 
¢, By differentiation of $(, y, 2, ¢, ¢,) = 0 with respect to x 
holding y constant there results 

Op O¢ 
On * Ba 
7 89 


p= 0. 
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By differentiation with respect to y holding x constant there 
results 


a ad 
y + ez dz 2 = 


By means of these two equations and ¢(2, y, 2, ¢; ¢) = 90, «, 
and ¢, can be eliminated. The result is an equation 


RGIS 2, ¢) = 9 
which is a partial differential equation of the first order. 


Exampie. Let a’+y’?4+2?+¢%+ ¢y = 0 be an equation 
between x, y and z, and two arbitrary constants ¢, andc,. By 
differentiation with respect to x holding y constant there results 


22 + ¢, + 2zp = 0. 
By differentiation with respect to y holding x constant there 
results 


2y + ¢, + 22q = 0. 


By elimination of ¢c, and ec, between the three equations there 
results 
ey — 2 Jeep 4 Qyeg = 0. 


This is a partial differential equation of the first order. 

(6) By the elimination of an arbitrary function. Suppose 
that w and v are functions of the variables 2, y and z, and that 
p(u, v) = 0 where $(u, v) is an arbitrary function of w and v. 

The differential of ¢(u, v) = 0 is 


od Od , 
Ln du An dv = 0. 
Now 
Ou ou 
du = aye + ay dz 
when y is constant, and 
du = ay as cds 


when 2 is constant, and similarly for v. 
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Therefore the partial derivatives of the equation ¢(u, v) = 0 
with respect to x and y, respectively, are 


Od[ Ou Ou oo Ov 
Ou dnt & | +5! Seta? | =9 
and ; ‘ 
p| ou ou Op¢[ Ov Ov 
4 | Se 1| +52 Spt oe2| = 9 
Eliminate a4 and 2 from these equations. 


; Ou ou Ov Ov Ov Ov Ou ou 
‘lant a? dy t dz! = \ dnt G2? (5; + a): 


When arranged in powers of p and q and the coefficients ex- 
pressed as determinants, the equation becomes 


|Ou ou (Ou Ou Ou Ou 
Oy Oz 2 da Ou Oy 
D = . 
au a0)" av dv)? ~ lav av 
| By Oz | Oz Ox Ox Oy 
This may be written in the form 
Pp + Qq =F () 
where 
ek eee ee 
Su dul Yow ul [ou day” 
Oy Oz Oz Ox Ox Oy 
ge ee eC 0% 5G 
Oy 02 Oz O02 Ox Oy 


This is a partial differential equation of the first order. There- 
fore from the equation $(u, v) = 0a partial linear differential 
equation of the first order can be formed which does not contain 
the arbitrary function $(¥u, v). 

ExampLeE. Suppose thatu=v2+y+4+zandyv=2#+4 y+ 2. 
Let o(u, v) =0 be an equation connecting w and v where 
$(u, v) is an arbitrary function of wu and v. é 
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By differentiation of ¢(u, v) = 0 with respect to x and with 
respect to y there result 


O O 
of 1 +p) + 52 (On + 2ep) = 0 


and 
oo od 
oy +9) + a (2y + 22q) = 9, 
respectively. By elimination of ae and = from these equations 
there results 
leans toe fake 
dS UB 
Zl Pee oreo 


or 
(y-2)p+ @—«)q=u—y. 


This is a partial linear differential equation of the first order 
which does not contain the arbitrary function 


oe+yt+4v+y +42). 


66. We have seen that a differential equation with two inde- 
pendent variables can be derived from an expression containing 
two arbitrary constants or from an expression containing an arbi- 
trary function of two independent functions of the variables. We 
see therefore that a differential equation with two independent 
variables may involve in its solutions, arbitrary constants or an 
arbitrary function of the variables. 


Definitions. A relation between the variables of a differential 
equation with two independent variables which includes two arbi- 
trary constants is called a complete integral of the equation. 

A relation between the variables of a differential equation with 
two independent variables which involves an arbitrary function 
of two independent functions of these variables is called a general 
integral of the equation. 

There is another class of solutions called singular integrals but 
these will not be considered here. 
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67. Consider the two equations u = ¢, and v = ¢, where u and 
v are functions of x, y and z, and ¢, and «¢, are arbitrary con- 


stants. By differentiation of w = ¢, and v = ¢,, there result 
Ou ou ou 
and y P 
v v 
respectively. . - 
. ) U 
Multiply (1) by Ap? (2) by a? and subtract. 
Ou Ou |\Ow Ou 
Oz Ox Oy Oz 
4 d= dy = 0. 
ov ov Ov Ov 
Goa lay bz 
Multiply (1) b ey (2) b ae nd subtract 
ultiply ( oe Lier act. 
Ou ou |Ou ou 
Ox Oy Oy Oz 
4 di | . dz = 0. 
Conte etl 
Ox oy Oy Oz 
dx dy - dz 
Ou Gul — [eu ul [au Ou 
Oy oz Oz Oz Ox Oy 
dv dv) Ja» ay| av ay 
Oy © Oz Ox Ox Oy 
Now ¢$(u, v) = 0 is a general integral of the equation 
FP Q R 7 
ae ou Su = ou Bu See Art. 65. 
Oy oz C25 C2 OL Cy 
go ov] ov av| fav av 
|Oy Oz Oz Ox Ox oy 
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Therefore $(u, v) = 0 is a general integral of the equation 
Pp + Qq = Ff if u =e, and v = «¢, are solutions of the equations 


68. From the investigations of Arts. 65 and 67 the following 
rule for finding a general integral of the linear partial differentia! 
equation Pp + Qq = # is determined. 

Solve the equations 


Suppose that w = ¢, and v = ¢, are two independent integrals of 
these equations. Then $(u, v) = 0 where (wu, v) is an arbi- 
trary function of wand v is a general integral of the equation 
Pp +Qq = B. 


Definition. The equations 


are called the subsidiary equations of Pp + Qq = R. They are 
also sometimes called Lagrange’s equations. 


69. As illustrations of the method of solution of a linear par- 
tial differential equation of the first order, consider the following 
examples, 


Examp.eE 1. Solve the equation 2p + ayq + y? = 0. 
Write the subsidiary equations 


Solve the equation 


Solve the equation 
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From 
Hi 
y = Cy, 


substitute the value of x, and the equation becomes 


A general integral of the original equation is therefore 
2 
& vy :) ZA 
EXAMPLE 2. Solve the equation (y — 2)p+ (z—a)q=2—y. 
Write the subsidiary equations 


dx dy dz 


EE ae a cama a 


From a familiar theorem of algebra, if 


then la + me + ne: 1b + md + nf =a: b,where 7, m and n are 
any multipliers whatsoever. Application of this theorem to the 
subsidiary equations gives 

dx + dy+dz=0, (1) 
when /= m= n, and 


ada + ydy + zdz = 0, (2) 


‘when 7 =2%,mM=Y,n=2. 

Solve equations (1) and (2). Therefore 7+ 44 2=¢, and 
“e+ty+t2=e, are solutions of equations (1) and (2), and 
therefore of the subsidiary equations. A general integral of the 
original equation is therefore (a +y+2,2+y? +2) =0. 
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EXERCISES 


Determine the partial differential equations of which the four 
following equations are complete solutions, ¢, and ¢, being arbi- 
trary constants. 


1. z= or + Gy. 2. 2 = 6,2 + ¢,y. 


on y 2 
3. 2=(#+64)4+4). 4, aitotaal. 
1 2 
Eliminate the arbitrary function from each of the four follow- 
ing equations. 


5. o(2+y—2z,¢7+y7—2)=0. 6. dba+ty+24 2) =0. 


to 2= OO(e-L ¥). eo. 2= te ye). 
Find a general integral of each of the following equations. 
9. xep — yzq = XY. 10. ep yg — 2 0. 
11. w’yp + yq = 22. 12. ep—yq =u —y. 


18. (' —2)p + @—a)q + (yt — 2) = 0. 
14) (2 — 3y)p + (Sa 42g = fy = oe 


ANSWERS 
1. ap +yq=%. 2. xp +yq=%. 
3. py = %. 4, wept yeq—2ta?= 0. 
o. (y—2)p+(2—2)q=y—x. 6. p—q=0. 
t. p—G= —2%. 8. yp —2#q = 


a i 
9. 6 (om logy + 9) =0. 10. +(5- 7 
1 ag ee 
U2. p(y, eo ye) 
13. d@+y+4e%4y742) =0. 
14. p(4e + 2y + 382, 27 4+ 47 + 2") = 0, 
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APPLICATIONS OF PARTIAL DIFFERENTIAL EQUATIONS. 
INTEGRATION IN SERIES 


APPLICATIONS OF PARTIAL DIFFERENTIAL EQUATIONS 


70. It is shown in the Analytical Theory of Heat that the 
change of temperature in any solid at a point (a, y, z) within 
the solid is given by the equation 

CU Os Wey Ou 

B=? | ae top toe ) 
where wu represents the temperature at the point and ¢ denotes 
time. 

In polar or spherical coordinates the equation becomes 


Ou Ou 
2 4 =, 
a(r S) 1 (sin 05% ) 1 du 


Ou 
é or sin 0 fale) sin? 6 O¢? ’ (2) 


ot 


and in cylindrical coordinates, 


Ou 4 E 1 Ou 1 Ou = | 


e 
Vad 


at" | OF TF or 1 A od? T OF ey 


Tf the solid is a rectangular plate so thin that the thickness 
need not be taken into account, equation (1) becomes 
Ou (ou. Ou 
EE LR ye Sees |" 4 
eae | Sart Se (4) 
If the solid is a wire of infinite extent so thin that the breadth 
or thickness need not be taken into account, equation (1) 
becomes ; 


= ee eee: 5 

at ~ ° Ox ©) 

In the case of a sphere when the temperature wu depends merely 

on the distance of the point from the center, equation (1), as 
97 
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can be seen from (2), reduces to 


O(ru) 20 ("U) 
ob. a ee (6) 


In the problem of permanent states of temperature, 0u/0t=0, 
and the equation becomes 
Cun, Cu 8% 
Ox’ ae oy Se o2? 
an equation known as Laplace’s Equation, and sometimes writ- 
ten V’u= 0. This equation also figures in the Theory of 
Potential. 

In polar or spherical coordinates, equation (7) becomes the 
right-hand member of equation (2) set equal to zero, and in 
cylindrical coordinates, it becomes the right-hand member of 
(3) set equal to zero. 

In the Theory of Acoustics, in considering for instance the 
transverse vibrations of a stretched elastic string, there occurs 
the equation 


0, (7) 


af = ° Ox” (8) 


and if the resistance of the air be taken into account, the 
equation 


ae t Fae = oa ee 


In the problem of the vibrations of a stretched elastic mem- 
brane, there occurs the equation 


Cee | Cee Ce 
ee © E or aa (10) 
which in cylindrical coordinates becomes 


Ove 2 Ore. 2 CL Oe. ke Oe 
Ce By ae arte O¢’ |" 


(11) 


71. As an ‘illustration of transformation of coordinates, con- 
sider the transformation of Laplace’s Equation in two dimen- 
sions, 
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Ou Ou 
02? =F oy ar 0, 
from rectangular to polar coordinates. 
The equations of transformation are x = r cos 6, y =rsin 6, 
Now w is a function of « and y and therefore of r and 6. 
Therefore, as seen in calculus, 


Ou Ou 
du u= Zs dr — a6 1. 


If y is held constant, this equation becomes, 


ou ou 
du = ap Og + a0 dO. 


Divide by Az, or what is the same, dx, and there results the 


equation 
Ou Oudr- Ou 0d 7 
Oa — Or dz 00 Ox’ a 
Similarly, 
Ou Ouodr Ou 00 
(2) 


Oy Or dy * 00 Oy 


Since « = rcos 6 and y = rsin 6, therefore r = Vx? + y’? and 


6 = tan y/x. 


or Ly x 
= — = 6080, 
oA ey 
Oy y =2=—sind 
Ve + y 
fale y Yad sin 0 
Ox ae RLY ciety oa 
and 
fala x c cos 6 
Cj a 
Now 


Ou Ouor ou OO Ou oer Ou sin 6 
Baie Cora Ob Of. OF Oy) er 


and 
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ou ofdéu es Ou sin 9 6 na Ou 9 Ou sin 6 sin 6 
Oat — Br| dr "BO -¥, Ore BO eA 
O*u j Ou snd Ousind Raye 
=~ | Or hl Orod r ni 06 r : 
Ou ip ola) 56 Cusind Oucosé|sin8 
late on + age er © 7 /eb) @ eee 
Similarly 
Ou Ou sin 6 O'u cos@ Ou cosé 8 
By | ar’ + ord. rf 8 «|S 


O*u 8 ee Cu cosO Ou sin 8 | cos 6 
T | B00 * Poet TOG ce 06. | 

, ou Oru > Cie LO Ores 

* Bat Cy or TF Or TF OO" 


Therefore the equation 


in polar becomes 


OUT SOu lee 
de tr ort Pog = 
72. A method of determining particular solutions of those of 


the above equations with constant coefficients is illustrated in the 
following example. 


Examp.e. Find particular solutions of the equation 


Oz Oe Oe 
z=" | zat ay 
Assume that there is a particular solution in the form 


z= e*tby+yt where a, B and y, are constants. 
Substitute in the equation. 


*, pewtbytyt — c(a 4 B?)exttByt 7, 
Now e***8y+v* cannot be zero for any values of a, y and t. 


Se y cas C(e ag ie. 
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Therefore z = e+8y+e’a*+8%t is a particular solution of the 
equation where a and £ are arbitrary constants. 

The above solution can be put into another form as follows: 

Let a = aj and B = fj where j = ¥—1. 


Then 
gg — eattByxctY a+ Bp»); 
Therefore 
z—sin (ot + By + ct va? + f°), (1) 
and 


2 = cos (az + By + ctVo? + 6”), (See Art. 5.) ~(2) 


are particular solutions of the equation. 
From these can be found particular solutions in the forms 


2 — sin ax sin By sin ct yo? + ’, 
z — sin ax sin By cos ct Va" + 8’, 
and six others. The determination of these six is left as an 
exercise to the student. 
73. Consider the equation 
Ou, 
O(r oe oO | sin gee 
or of 1 CO | 0 (1) 
or sin 6 fale) ere 


which is Laplace’s Equation in spherical coordinates where w is 


independent of ¢. 
Let u = r"P, where P is a function of 6 alone, and m is a 
positive integer. On substitution there results the equation 


i a( sin Oe ) 
m(m + 1)r"P + = = 0: 


né 00 
Change the independent variable from 6 to x where x = cos 0. 
he ad dP 
: BA ie wet BO gta! = 
at i pe 2x 7 + mm + 1)P= 0. (2) 


The solutions of equation (1) are known when P is determined 
from equation (2). Equation (2), not only when m is a posi- 


\ 
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tive integer but for all values of m, is called Legendre’s Equa- 
tion. Its solutions are discussed in Arts. 76, 85 and 86. 


74. To find particular solutions of the equation 
Cree ee: lee 
ge = 9 [Sats orf © 
which is equation (11), Art. 70, when z is independent of ¢, 
let 2 = R- T where RF is a function of 7 alone and T is a fune- 
tion of t alone. Substitute in the equation. 


eT ES a) 


yt de" r dr 


Sie 


or 
(2) 


The right-hand member of (2) does not involve t. Therefore 
the left-hand member does not. The left-hand member does not 
involve 7, Therefore the right-hand member does not. There- 
fore each member is constant. Call the constant — yp’. 


Lad road E 1 i | 


Di dit Paleo ae 


fife 
cee +peT = 0, (3) 
and 
abla AN Ales 
det, gp tht =. (4) 


A particular solution of (1) is therefore R- T where T is de- 
termined by equation (8) and R by (4). 

Particular solutions of equation (8) are J’ = cos pet and 
T =sinpct. (See Art. 31.) 

To solve equation (4), let r= «/p and substitute in the 
equation. 


Cen leroy 
. 2 . dv + v dz + R = 0. x (5) 
Equation (5) is a special case of the more general equation 


dl”, dt 
fates + (2 —n*)y = 0, 
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known as Bessel’s Equation. Its solutions are considered in 
Arts. 79 to 84 inclusive. 


INTEGRATION IN SERIES 


75. It will be noted that as yet in this book no equations with 
variable coefficients, of higher order than the first, have been 
considered except a few very special cases discussed in Chaps. V 
and VI. The remainder of this chapter is devoted to a discus- 
sion of linear differential equations of the second order with coef- 
ficients rational integral functions of «, and second member zero. 
To such a set belong Legendre’s and Bessel’s Equations men- 
tioned above. 

Not all differential equations, not even all in the comparatively 
simple form of linear differential equations of the first order, are 
capable of solution in finite form. When solutions cannot be 
found in finite form, recourse is had to integration in series. In 
the set about to be considered, some equations have solutions in 
finite form and some have not. 

We shall attempt here to find solutions only in the form of 
infinite, convergent, power series. 

If an equation be capable of solution in finite form, this form 
is found when a solution is attempted in the form of a power series. 
For instance, in exercise 11, page 120, the solution found as if it 
were made up of infinite series is in reality in finite form. 

Sometimes the series that make up the solution of an equation 
may be recognized as those of familiar functions. In such cases, 
the solution can be written in terms of those functions. For 
instance, in the answer given on page 122 for exercise 12, page 
120, if A be taken equal to 2 and B to 1, and the two particular 
solutions be added, there results the series 

(2x)' | (20)* (20)" 


which is w# %e”. If the oe solution be subtracted from the 
first, there results the series which is 2 *e**, The general solu- 
tion is therefore 


uy 
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Y= OH 7e* + Or 6. 


76. Let us attempt to find a solution of Legendre’s Equation 
2 Uy dy 
(1 — #7) 54 — Aas + mm + 1)y = 0, 


in the form of a power series in 2. 
At first, assume that there is a power series 


Pn Di a Pecos et oa eo pee aa 


where g,, 9, °°, « are constants, which will formally, i. e., 
without regard to whether the series converges or not, satisfy the 
equation. It is no restriction to assume, as we shall, that g,=-0, 
because, if there is any solution at all, one at least of the g’s is 
not zero, and we assume that the series begins with the term con- 
taining the first g which does not vanish. 


Since 
y= Zoe 
. dy —e ak y— rer 
ae 
and 


ce 


=> (x0) (eel gee aes 
Substitute in the equation. 


< a [1l— a) yt val ga 
— Qaln + v)g,arte— 4. mim 4 1)g9,0*t"] = u 
or 
oD [ke + vy) («+ y— T)gar 
Zz { (k +} v) Ce ee ee 1)}g,aet"] ome (1) 


v=o 


ye yo Oe 


v=0 
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is to satisfy the equation, the coefficients of each power of x in 
(1) must be zero. Therefore there results the following series 
of equations : 


K(k—1)g,=0, (x+1)«g, = 0, 
(+ G+ 1g, — {x(« + 1) — mem + 1)}9, =0, 
(x + 3)(« + 2)g, — {(« + I(x + 2) — mim + 1)}g, = 0, 
(2) 
Gee arn 4-2 — 19, 
= {(« + 2r — 2)(« + 2r —1) —m(m+1)}Q,_, = 0. 


From the first of these equations, since g, + 0, therefore 
x=0, ork =1. At first, take « = 1. 
Substitute in equation (2) and calculate the g’s in succession. 


Soi | 2 
9, = 9, eee tea doy 9, = 9, 
el 9 2 4 
T= o un ie —— ae ah, Gang = 0, 
ety 
sm 1)(m—38)---(m—2r+1)(m+2)(m+A4)-+- (m+ 2r) 
j2r 41 Jo» 


(m—1)(m +42) 


eee | @ 


where g, is arbitrary, and g,, has the value given above, formally 
satisfies the equation. Since this series is convergent, (3) is a 
particular solution of the equation. 

Next, take x = 0. 


Substitute in equations (2). g, is arbitrary. Call it zero. 
8 t 


106 SHORT COURSE ON DIFFERENTIAL EQUATIONS 


m(m + 1 m(m—2)(m+1)(m+3) 
p= oy = 0 = |4 0 
pie F549) 
, m(m—2)-+-(m—2r4+2)(m41)---(m42r—1) 
I= (—1) ( ) |2r Iw 


m(m-+1) m(m—2)(m+1)(m+38) 
y= 9/1 - 2 a? E oe 
4+ 95, gr ue. | (4) 
Io 
where g, is arbitrary, and g,, has the value given above, formally 
satisfies the equation. Since the series is convergent, (4) is a 
particular solution of the equation. 
If solution (3) be denoted by y, and solution (4) by y,, the 


general solution of the equation is y = Ay, + By, where A and B 
are arbitrary constants. (See Art. 11.) 


77. The general form of a linear differential equation of the 
second order with right hand member zero is 


dl? d ’ 
go(2) Fa +.4(2) 5! +4,(2)-y = 0. (1) 


- It will be assumed here that q,(x), 9,(@), q,(#) are rational 
integral functions of x. 
If a solution is-to be found in the form of a power series in 
x — a, it will be convenient to write the equation in the form 


(@- ope) 21+ @—ap@ L+7,@)-¥=0, @) 


where p,(), p,(#), p,(#) are rational integral functions of 2. 
The equation can be written in this form in an unlimited number 
of ways by multiplying it through by a suitable power of x — a 
and a rational fraction neither the numerator nor denominator of 
which contains « — a. 
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Definition. The point a is a regular point of equation (2) if 
Po) + 0. 


Without at first making any assumption with regard to the 
point a, substitute 
eee 2) 


equation is formally satisfied 


v= 


D3, 


in equation (2) and attempt to determine the g’s so that the 


[Oty Drv) +) 2) +) [ov @— ye 


=0. (3) 
Call the expression in square brackets f(x, x + v). 

Develop f(x, x + vy) into a power series in « — a by Taylor’s 
Theorem. 


fa K+) =flan+v) + fila «ty = Ti = 4 


a Cae eee 
nv 
Substitute this development in (3), equate each power of 


x — a to zero and there results the following series of equations 
Iof(4, x) = 9, 
nia, K+ 1) + 9 f' (a x) = 0 
, : i" (a, K 
9.F(4, «+ 2) + HF © +1) + 9% = FUG) 


Gila, En) +9, 4) (4, Kn —1) 42 + 4+ fe oe 


ae 
f(a ntl) _ 4 
m+1 ae 


Invid (4 Kant 1) gf’ (a, Kalai inte = to 


Now 9, + 0. 


Therefore f(a, «) = 0. 
And 


S(4 «) = K(k — 1)p,(@) + «p,(4) + p.(@) 
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5 K(K ou 1)p)(@) =F Kp, (a) + p,(4) = 0). 
From this equation can be determined the value or values of x 
which are to be used in the subsequent equations (4). If p,(a) 
is not zero, the equation is of the second degree. If p,(a) is 
zero, the equation is of lower degree than the second. 
The necessary and sufficient condition that the equation 


K(k oF 1)p,(«) = Kp, (a) + p,(@) =.0 


is of the second degree is therefore that the point a be a regular 
point of the differential equation. 


Definition, The equation x(x—1)p,(a)+«p,(@) + p,(a) = 0 
is called the indicial equation of the differential equation (2). 

If the point a is regular, the indicial equation gives two values 
of x, say x’ and «’’, and from equations (4), for either value of 
x, the values of g,, g,, +--+, may be computed, in general, in 
terms of ,. 

Therefore in general there are two series in ascending powers 
of « — a, namely, 


y= = 9, (@ —a)'*”, 


and 


y= Lgp(e—ay'”, 


where g, is arbitrary in either series, which will formally satisfy 
equation (2). 

78. The following theorems with regard to the solutions of the 
differential equation 


(a — a)'p(2)54 + (@ —a)p,(o) Y + p(x) -y = 0 


in a power series in # — a have been established. The proofs 
are too long to be given in this book. For a discussion of these 
theorems the student is referred to a pamphlet entitled ‘‘ Regu- 
lar Points of Linear Differential Equations of the Second Order’? 
by Professor Maxime Bécher, published by Harvard University. 
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Theorem I. If ais a regular point of the differential equa- 
tion, and the difference of the roots of the indicial equation is not 
zero or a positive integer, two solutions in the form of a power 
series in x — a, viz., 


y= 


p= Ener, 
and “ 


y= Fae aym 


where x’ and x” are the roots of the indicial equation, exist, and 
these series are convergent. In each of these series g, is 
arbitrary. 

In this case, if y, denotes one of the series and y, the other, 
the general solution of the equation is y = Ay, + By, where A 
and B are arbitrary constants. 

A case to which this theorem applies is Bessel’s equation when 
n is not zero nor an integer, discussed in Art. 80. 

Note. By the difference of the roots of the indicial equation 
being a positive integer is meant that the greater minus the less 
is a positive integer. 

Theorem II. If a is a regular point of the equation and the 
difference of the roots of the indicial equation is a positive inte- 
ger n, the necessary and sufficient condition that two solutions 
of the form under Theorem I exist is that 

(n) ” 
GRRE) 0, 


jn 


Gat’ (a, KT to — 1) oes 4g, 


(see equations (4), Art. 77), where «” is the smaller of the 
roots, and when this condition is fulfilled, the series are conver- 
gent. 

In this case the series corresponding to the larger value of x, 
say x’, can be found as before. In the series corresponding to 
«’, g, and g, are arbitrary, but if g, be chosen zero a particular 
solution in terms of g,is found. Then if y, denotes the first 
series, and y, the second, the general solution of the equation is 
y = Ay, + By, where A and B are arbitrary constants. 
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A case to which this theorem applies is Legendre’s Equation 
discussed in Art. 76. 


Theorem III. Ifa is a regular point of the equation, and the 
difference of the roots of the indicial equation is either zero, or a 
positive integer n where 


" (n) a, wi 
Geen Cae $n a1 po tt 0, 


two solutions are found, one being 

yon 
Y= Ve — ayer”, 
: v=0 
the other being 


y = log (« — a) Loo = ee LHe — ayer”, 


where «’ is the larger root of the indicial equation and x” the 
smaller, and these series are convergent. In the first of these 
series g, is arbitrary. In the second, g, is arbitrary and g, is 
determined in terms of 9). 

If y, denotes the first series and y, the last term of the second, 
the general solution of the equation is 


y = [A + Blog (« — a)] y, + By, 


where A and B are arbitrary constants. 
A ease to which this theorem applies is Bessel’s Equation 
when n = 0, or an integer, discussed in Arts. 82 and 83. 


Theorem IV. If the point a is not a regular point of the 
equation there are not two solutions of the equation in any of 
the forms under Theorems I and III, and if any series in one of 
these forms is found it is usually not convergent. 

Cases to which this theorem would apply will not be considered 
in this book. 

BESSEL’S EQUATION 
79. We shall now consider the solutions of the equation 


d*y d 
ab eae + (2 — n*)y = 0, 


in the form of a power series in 2. 
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The equation as it stands is in the form (2) of Art. 77, where 

pl«#) =1, pv) =1, P(e) =v — Nn. 

Since p,(x) cannot be zero for any value of z, all points of 
this equation are regular. Therefore the solutions of the equa- 
tion for all values of x, and in particular when « = 0, will come 
under one or other of the forms mentioned in the first. three 
theorems of Art. 78. 

Substitute 


v=n 


] a a Opp rre 


in the equation. 


Pa: S [«+rvy(e+v—1) 4+ +) 4+ (—2’)] gat" =0, 
or 
be [«+¥)? +0 —n']g,axt” = 0. 

The equations for the determination of the g’s are therefore : 
(x? — n*)g, = 9, 
[(«x +1)? — n’]g, = 0, 
[(« + 2)’ — n*]9, + 9 = 9, 

Eas — I. = °, (1) 


ice + or — 1) — alms re 0, 
[(k Se 1 “ez VI ap Io, = 0, 


Since g, + 0, from the first equation there results x = =n. 
The difference of the roots of the indicial equation is therefore 


+2n. If n =0, this difference is zero. Ifn = see » where p 
is an odd integer, or if 7 is an integer, this difference is a posi- 
tive integer. If is neither zero nor +5 nor an integer, the dif- 


ference is neither zero nor a positive integer. 


2 


80. At first assume 7 neither zero nor = 9 nor an integer. 
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This is the case covered by Theorem I. There are therefore 
two solutions 


ee au war 
and 


v= 


y= zu ie ene 


To determine the g, substitute x = n in equations (1). 


Io 


>9,= 0, 92 = — (In 4 2)? 9, = 9, 


Geum I 
Is = — Ain +4) 2-402n + 2)Cn £4)’ 


Goy—1 ie 0, 


Io 
I = (— WY 5g 


- Qr(2n + 2)(2n 4+ 4)--- (2n + 2r) : 


ae tf ache oe ae 
: 2( 2% +2) 


4 


5 4Gn £2)n 4) oe | @ 


where g, is arbitrary and g,, has the value given above is a par- 
ticular solution of the equation. 
Similarly, on substituting « = — in equations (1) there 
results 
aid 


Y= oe” E +50 2) 


at pt es arene 
+9. 40n= 2) Ona 4) segs | 


where g, is arbitrary and 


ae Io 
wi 224s (Qn — 2) Ome) Sana 


is a particular solution of the equation. 
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If y, denotes the first series and y, the second, the general solu- 
tion of the equation is y = Ay, + By, where A and B are arbi- 
trary constants, 


81. Next, assume n = 2 Assume, for definiteness, that p 


OM 
is positive. In this case the difference of the roots of the indicial 
equation is a positive integer, viz. p. From an examination of 


equations (1) it is seen when n = - that both g, and g, are 


arbitrary. Choose g, = 0, and there results the same equation 
as (3) of the preceding article when +5 is substituted for n. 


Therefore in-this case there are two particular solutions of the 
equation which are the same as the solutions in the case of the 


preceding article when ? is substituted for n. 


2 
82. Next, assume m an integer. Since n appears only in the 
form of a square in the differential equation, it is sufficient to 
suppose it a positive integer. 
In this case the difference of the roots of the indicial equation 


is the positive integer 2n. 
For the root x’ = n, the series is the same as (2) of Art. 80. 


For the root «’” = — n, the equation 


[(- Nest 2n)? a 1”) 9Ion + Yn-2 = 0 
is such that the coefficient of g,, is zero. Therefore, since 
Gon» + 9, this case comes under that mentioned in Theorem ITI, 


Art. 78. 
To get a solution corresponding to «’’, let 


y = loge D gyar’ty 4D Gyan”t”, 
v=0 v=0 


For the purpose of determining the coefficients g,, write the 
series in the form 


y=n 


ee au (Gi NE SIG oe 


where 91, = 9m = '°* =9.i=°%. 
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Substitute in the equation. 

ne > [(R” ey)? a? = 90 ge, le eae” 
v=0 


+ {2K 4-7) Gyan [K+ v)? +0? 7] 9, ort =0. (1) 
Now in the coefficients of log x- x*"+’, 9_.., Yoniyy °° *s Y--» 
are zero, and the remaining ones are the same as the left hand 
members of equations (1), Art. 79, with «’’ + 2n substituted for 
xk. Now «”’ 4+ 2n =x’, and equations (1), Art. 79, hold for 
« =x’. Therefore all the coefficients of log x - x*”+” vanish. 
From the second set of terms in (1) are found the equations 
from which to determine g,. These equations are: 


Comer Dre ie 
[(«” 4 1)* — n4]9, = 9, 
[(«” + 2)* — n*] 9, + J = 9, 


[Cx + Qn)” — 27} Gan + Gina + 2(K" + 2n) gq = 0, 
[(«" + 2m + 1)? — 2) Gansr + 2K" + 2n + 1) g, = 9, 
[CK + Qn + 2)? — 27) Oynse + Jan + 2(K"” 4+ 2n + 2)9, = 0, 


Since in these equations, «’’ — — n, therefore g,, an arbitrary 
constant, satisfies the first equation. Also, 
Re) gee 
oS OM 58s Cea 
9. = 0 Gg = Jo Aah 
den I= Aan = Oe 


Io 
-27r(2n—2) (2n—4)- - -(2n—2r) 
Since the coefficient of g,, is zero, this equation introduces no 
new g. The equation, however, gives a means of determining the 
hitherto undetermined constant g,, 


921=9; r<oN, Ix = 5 ae 


Ces es 2ng,. 
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Also 
g = Jo 
S— 5) aes | Winer 1 1 
Jp = — 2" |n|n—1y,, 


Io 


A SRN ie ee, a ea 
n= b= Fim ay °° Sma = TE 
where 

Jo = — 2" |n|n—1y,. 
Choose g,, = 0. Therefore 


Daas = 0, 


“Qa = 2)}" 


4 _ _ (m42)9,_ (+2)q 1 [; 
ae wean ee Op) Deon?) | 2 


Iontara = 0, 


1 


il 
tar {ver 


a= (-1)""5 Re Gee 2r(2n42)(2n+4) tae 


(Qnp2r) 7 


aul itt il 
Ko are oe) 


ees oe = Ki+y __ 92n-1 = te TO ON 
F NRE 2 [nn —1 gg E + 5(in — 2) 


i x 


ay? 


+9-4(Qn—2)(Qn—4)* °° (ceca 


ge il 1 
Ber) 


9 Sie on + 2} ocr Jontare ie eae 


where g, is arbitrary and g,,,,, has the value given above, is a 


particular solution of the equation 


If the first solution be denoted by y, and all terms not involv- 


ing log w in the second by y,, the general solution 


of the equation 


is y = (A+ Blog x)y, + By, where A and B are arbitrary 


constants. 
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83. Next, assume n = 0. Im this case the difference of the 
roots of the indicial equation is zero. This case is covered by 
Theorem ITI. 

The two infinite series when n = 0 can be found from equa- 
tions derived as in the preceding case. They can also be found 
by letting n be zero in the results in that case. Two particular 
solutions are 


a? at qr 
Pr 1 — pe + oF Bot) Fn 
and 
x? at x feat | 
as Comey — fi DN ae eee Eo 2 
Ya= 9 log *+ 55 +a get +( ty FP. & r 


84. As will appear in applications to physical problems, when 
n is a positive integer it is convenient to take, not y, and y,, but 
the quotients of these by 2"|n where g, is unity. These special 
solutions are written J,(2) and W(x) so that 


qr x? at 
JA) = 9575 [1 mney +7 aC Cree | 
T=0 (= Lee 


= 2X te |r ee 
and 
Ors Nee Da nt? 
W(x) = J,(x) log x — a" n—1le” + aie 
ore ”— 8 gtr gn? 
ar |2 +b 


r= (= Rigas”: nt+2r rr 1 all 
zr pe af Ortr |p rln+r 2 (s+ sa) 


LEGENDRE’S EQUATION 


85. Returning now to the equation 


(1 aoe Se 20 il +m(m+1)y = 0, 


considered in Art. 76, we see that it can be transformed into the 
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form (2) ‘of Art. 77 by multiplying through by 2%. The trans- 
formed equation is 


d? 
Celia ) “ia =o be + m(m + 1)a’y = 0, 


where p,(z) =1—2', p,(#) = — 227, p,(#) = m(m + 1)z’. 
Since p,(#) = 0 when x = — 1 and 2 =1, the points — 1 and 
1 are not regular points of the equation. All other points are 
regular. 

The indicial equation when # = 0 is 


K(k —1) = 0. 


This equation has the roots 0 and 1. The differential equa- 
tion in this case comes under the case mentioned in Theorem IT. 
This case was already discussed in Art. 76. 


86. It is convenient when m is a positive integer as it was in 
the illustration of Art. 73 to have a solution of Legendre’s Equa- 
tion as a series in descending powers of x. In this case we shall 
take the equation, as in Art. 73, in the form 


or dP 
Ch) ie 22-7, +m(m+1)P=0. 


Let 
D Se iO Pi ae 
v=0 


and substitute in the equation. 


e a (n —~v)(n —v —1)g,9°-"-2 
v=0 SiG 1) na yb) min 41) ga = 0. 
a inGa Ll) mcm 1) gy = 0, 
[(n —1)n — m(m + 1)]9, = 0, 
n(n —1)g, —[(r — 2)(n —1) — mim + 1)]9, = 0, 


(n — 2r — 2)(n — 2r — 3)g,,_, 
— [(@ — 2r)(m — 2r + 1) —m(m + 1)]9,, = 9. 
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From the first equation, since g, + 0, 
n(n 1) — mn y= U; 
“.n=m or n= —m—t1., 
“ft first take n = m. 
(m —1)m 
= 974=9 %=- 2(2m — 1) 9 9s=9, *** Drea = 9, 


1 ue Saha Oo Cae — -(m — on 
at A°2 Om = Wr ly: “One iy 2 


A series in descending powers of « which satisfies the faaaae 
is therefore 


* (m — 1)m pe 
aoe j1- 2(2m = ete 


(m — 8)(m — 2)(an — 1)m - Gs a ae 
Bo AO (Ow 2 5) (ee ee | © 


~nere g, is arbitrary and g,, has the value given above. 


By taking n = — m —1, there results the solution 
Lace (m + 2)(7 +1) = 
a ae E + 9am + 8) 
(m + 4)(m + 38) (m + 2)(m + 1) gta ss Seg (95 
4-2(2m + 38)(2m + 5) TG, 


where g, is arbitrary and 


Rs aie LE ad -(m +1) 
— Or.--2(2m + 8) - oni + 2h) Jo: 

When i is a positive integer, solutions (3) and (4) of Art. 
76 is a finite series according as m is even or odd, and in either 
case, equation (1) above is finite differing from (3) or (4) of 
Art. 76 only by a constant factor. 


87. If series (1) of Art. 86 be multiplied by 
Cn) m3 ard ne |2m, 


[m 2"(/m)* 


the resulting integral is 
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called the Legendrian Coefficient of the mth order, and is de- 
noted by P (2). 
The successive values of P(x) are readily found to be 


Pat Playas, Px). $0 — 4, 


EXERCISES 


1. Find the remaining six particular solutions of the equation 
considered in Art. 72. 
2. Find two particular solutions of the equation 


oN 


C= on! 
3. Find two particular solutions of the equation 


Oru) eo Or (ru) | 


4, Find four particular solutions of the equation 
Ou Ou 107 
ot \ Bat + By 


5. Find four particular solutions of the equation 


6. Find four particular solutions of the equation 


ory Oy OY 
an + hae = © au 
7. Show that equation (1) of Art. 70, in rectangular coor- 
dinates, becomes equation (2) when transformed to polar or 
spherical coordinates. The equations of transformation are 
z—rcos@sin¢g, y=rsin@sin dg, z=7 cos ¢. 


ou Oudr dOucd du Od 


Oe imilarl 
pA ce eo ee ey 


Suggestion. 
for y and z. 
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Find the general solutions of the following equations: 


dy dy 
2 St! 2 = 
8. 2 Aes jg tp L+P)y = 
a dy 
9 9a 5 + (92 — 2") = Cp 0: 


10. 4 Pod, dob (4 ay = 0. 


11: ood _ eS! + oy = 0. 


dy d? 


19, 54 4 40 V+ (2—42")y = 0. 13. i — ay = 0. 


14. (@— 2% Loe t+ 2y = 0. 


a 
15. « ae + Qy = 0. 16. 2f ey 0. 


17. Show that, 
Tiiy) = —J,(y), 
2, (y) = J,_4(Y) aa J Ws 


Qn 
yy J,,.(y) = Jay) ae Tri). 


ANSWERS 

1. z= sin aw cos By sin ct la? + 6, 
2 — sin ax cos By cos ct Va? + 8’, 

2 = cosaasin By sin ct Vo? + 8’, 

z = cos ax sin By cos ct yo? + 8’, 

z = cos ax cos By sin ct Vo? + 8’, 

z = cos ax cos By cos ct Va? +. f°. 


2. uw=e—%* cosax, U = et sin ax, 
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pf ©. if ° 
3. uU= 3 es cos ar, U = = e~a* sin ar. 


4, u = eH +B% cos ax cos By, U = e~*(a°+8) cos az sin By, 


U = e+ BY sin ax sin By, w= e~*(2?+ 87 sin ax cos By. 


5. = COS a& Cos cat = sin ax sin cat 
? ? 
= sin a& COs ca = COS ax sin cat. 
n O i, =UCO} t 
A =e Ca = é sin ax sin Ca —K 
6 ee Beye iss t "OQ? hes 


y = €™ cosax cost Veu—F, y = &cosarsint Ver? — i. 


1 ae 
8. yods|t— oat tras" 


: 5 od CL BE | tae 
oes Uris PPR [+2 eee 


0 eer : Of +++ 4 9,0" 
Saletan POA Gt LL Sor : 
where 
ie : 
In =) 2-4. 1Qr-5-+9-+- (Qr41) 
and 
— 1)” 1 
9x, = (—1) Dede Dy Bie. (Apo): 
oe 


7 4.7-10 
eg lt+55 e+ R56" +355.16.337 tc t9e + 


2:5, 2-5-8 ae 
+ Ba- t+ et gra + 36.917 an 


where 


7: +4 3r) a. 2+: @r—1) 
ee, rt Br I) 


10. y = Aci E teat tc tet | 


at. 


Blt par rasta + Gat], 
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where 


1 <I 1 
Is = F210 (IP gor) oe eee 


Wh, y= Ae Be. 


= 1 as 
12, y = Ax J. + fe + gage + Sart” +: =| 


2 


4 se 
+ Ba” E $4 tye te the +], 


4" = 4” 


= G0 rar) oo = a eee 


8% y= del 1+ se is 


1 
1g? SEA, Fae sou 
where 
: nd 4g. : 
#19749 ---3(87 27) a ju = 3) age 
14.0 y = Al — 40 oe) Ae Be [1 4a], 
15. y= [A+ Blogz] 


bo 


2 2 2 A . 
x 1- +77? — 7F.o8. ct ae 


Ors 98 De 
+B| he — pig ®C +2) + mong? +343) 


vhs re 
—jqop gr pt +2+34+4) tet Ge te, 
where 
QF grt r=r | 
9, = (— 1)! ine and 9. = (1) ae 


(eth 
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2 3 aS 
16. y=(A+Blogz)}1— gv + a pt: Ba Pe ae le 


2 


we = 2 2 
+3) 8 ope $8) te peed th4 


a 
where * 
2a Cee GE al 
ee FL (Ory? aE Miers, Big! UO; me 


FISHER AND SCHWATT’S 


Series of Text-Books on Mathematics 


Rudiments of Algebra...........$ .60 School Algebra...... ssceeeeeee $3.00 
Secondary Algebra.........-....- $08 Elements of Algebra.........-.. 1.40 
Complete Secondary Algebra.. 1.35 Higher Algebra..................- 4.50 
Quadratics and Beyond......... .90 ‘Text Book of Algebra, Part 1. 3.40 


FHE Fisher and Schwatt algebras have achieved a marked success 

in the short time that they have been before the public. The 

‘“‘ Higher Algebra”? is used as a regular text-book at Harvard 

=} University, and the elementary books are established in sec- 
ondary schools of corresponding rank. 

These text-books may be resolved into two series, the more popular of 
which consists of the ‘‘ Rudiments”? (for grammar schools) and the ‘‘ Sec- 
ondary’”’ (for high schools). The ‘‘Complete Secondary’? includes the 
material of the ‘‘Secondary’’ and also chapters on Continued Fractions, 
Summation of Series, the Exponential and Logarithmic Series, Determin- 
ants, and Theory of Equations. Thus, while the ‘‘Secondary”’ covers 
the requirements for admission to colleges, the ‘‘Complete Secondary”’ 
covers the requirements for admission to any scientific school, and is also 
sufficiently full for the ordinary work of the first yearin college. ‘‘ Quad- 
ratics and Beyond”’ consists of the second half of the ‘‘ Complete Sec- 
ondary,’’ bound separately for the convenience of advanced and reviewing 
classes in secondary schools, and for college freshman work. 

The second series constitutes a more difficult course and places more 
emphasis on the theory of mathematics. In this series the ‘School 
Algebra”’ corresponds roughly to the ‘‘Secondary,’’ and the ‘‘ Elements”’ 
to the “ Complete Secondary.’’ 

The ‘‘ Higher Algebra’’ isa book for college courses. The ‘‘Textbook 
of Algebra’’ is a high school book containing an unusually large number 
of graded exercises; it is a valuable mine of problems bearing on the 
regular high school work. 


THE MACMILLAN COMPANY 
64-66 FIFTH AVE., NEW YORK 


BOSTON CHICAGO SAN FRANCISCO ATLANTA 


The Elements of the Differential 
and Integral Calculus 


By DONALD FRANCIS CAMPBELL 
Professor of Mathematics in Armour Institute of Technology 


1zmo. Cloth. x+ 364 pages. $1.90 net 


HIS book is designed especially to introduce the student 
of engineering to the mathematics upon which his future 
ORAS) work will be based, but in spite of the emphasis on the 
2 practical side of the subject the needs of classes in clas- 
sical colleges and universities have been neglected in no way. 

To meet the needs of colleges where the study of the calculus 
is taken up in the first year of the course, Professor Campbell has 
presented a more detailed discussion in the opening of both the 
differential and integral parts of his work than is usually given 
in text-books. Thereafter, the subject is developed by the use of 
practical problems which are sure to arise in engineering work. 
Thus all subjects only remotely connected with engineering have 
been omitted, while in addition, a few elementary chapters in 
mechanics have been supplemented. This presentation of mate- 
rial, without encumbering the book, affords a short introduction 
to Mechanics and Differential Equations as well as a view of the 
principles of Attraction, Centers of Gravity, and, to a certain ex- 
tent, the Moments of Inertia, from the mechanical rather than 
from the purely mathematical side. 

The part of the book which differs most widely from other 
text-books is that dealing with the integral calculus. A full ex- 
planation is given of each step in the formation of each summation 
and integral. In addition, in order to enable the student to grasp 
more fully the details of the subject, the author has introduced a 
large number of practical questions which are found in actual ex- 
perience to produce the desired result better than the theoretical 
propositions introduced into the older treatises, 
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